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1.0 


SUMMARY 


This study was performed under NASA Contract NAS1-14031 to 
define a configuration for the B-52E aeroelastic wind tunnel 
model that would produce high-frequency symmetric and anti- 
symmetric flutter modes with violent onset, synthesize a 
flutter suppression system (FSS) capable of stabilizing 
these modes and implement the FSS using digital computers. 

The system was then tested in the Transonic Dynamics Tunnel 
at NASA Langley Research Center. 

For the past ten years the Boeing Military Airplane Company 
has assisted NASA Langley in demonstrating the feasibility of 
active control systems designed to augment or suppress low 
frequency structural modes. During the B-52 CCV program a 
1/30 scale B-52 aeroelastic wind tunnel model was used to 
predict and verify the performance of the ride control and 
flutter suppression systems. The results of the wind tunnel 
and flight tests showed good correlation to analytical results, 
verifying the usefulness of this methodology. 

The flutter mode on the B-52 CCV program was a symmetric, 
low-frequency (2.4 Hz) mode with mild onset. For the current 
study, the B-52E aeroelastic model was modified to produce 
symmetric and antisymmetric flutter modes with violent onset 
at higher frequencies (13 Hz to 25 Hz model frequencies; 

2.4 to 4.6 Hz airplane frequencies). A flutter suppression 
system was then synthesized to extend the flutter dynamic 
pressure of the modified model at least 44 percent. The 
resulting FSS feedback filters were mechanized using digital 
computers with three channel redundancy to provide fail -operate 
capability. 
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2.0 


INTRODUCTION 


This document Is the final report of Contract NASI -14031 entitled 
"Stability Augmentation System for Aeroelastic Wind Tunnel Models. 

The primary result of this program was a flutter suppression system 
Implemented digitally which was capable of stabilizing the flutter 
modes of the modified B-52E aeroelastic wind tunnel model. This 
system was successfully tested In the Transonic Dynamics Tunnel at 
NASA Langley Research Center. 

The work began with the definition of the structural changes required 
for the existing l/30-s1ze, full-span, cable-mounted, free-flying 
model of the B-52 CCV airplane to obtain the flutter characteristics 
required for this study. (Descriptions of the original B-52 model 
and previous wind tunnel tests results are presented In references!, 
2, 3, and 4). For the modified model, flutter was to consist of two 
modes, symmetric and antlsynmetrlc, that exhibited violent onset and 
had frequencies In the range of 13 Hz to 25 Hz (model frequency). 

To allow testing above the flutter speed without exceeding the design 
limit of the model the flutter modes were required to occur below 
3831 N/m2 (80 psf). A configuration which satisfied these require- 
ments was defined that had wing ballast replacing the engine nacelles 
and external fuel tanks and was sting mounted at the wing attach 
point as Illustrated on Figure 2-1. 



FIGURE 2-1 - PLAN VIEW OF B-52E AEROELASTIC MODEL 
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The flutter suppression system was synthesized to produce a 44 
percent increase in flutter dynamic pressure with ±6 dB gain 
margin and ±0.7854 rad (45 degrees) phase margin at speeds below 
the flutter velocity. Due to the similarity of the symmetric 
and antisymmetric flutter modes the final symmetric and anti- 
symmetric feedback filters were identical, as shown on Figure 
2 - 2 . The final system was predicted to extend the flutter • 
dynamic pressure more than the required 44 percent as illus- 
trated on Figure 2-3. 

The FSS was implemented using digital computers in a three 
channel, redundant arrangement. The filters were transformed 
from analoq form into difference equations using the bi -linear 
transform (Tustin's method) and implemented in parallel form 
using scaled integer arithmetic operations. The system per- 
formed failure detection using a circular comparison technique 
where each computer compared the voter output to another com- 
puter's output. The final form of the FSS is shown on Figure 
2-4. 

The wind tunnel test of the model equipped with the FSS was 
conducted in 95 percent freon with a density of 2.58 kg/m'^ 

(0.005 Ib-sec^/ft^) at 6400 m (21 000 feet) equivalent air- 
plane altitude. Although the flutter dynamic pressure was 
considerably above the predicted value, 3926 N/nr (82 psf) 
instead of 2873 N/m^ (60 psf), the FSS performed as predicted, 
stabilizing the flutter modes up to 4884 N/m^ (102 psf). In 
addition, the FSS was evaluated with various induced failures 
and degradations. These tests proved the fail -operate capability 
of the FSS and the ability of the system to reduce the effects 
of a channel degradation when another channel had failed. Post- 
test analysis revealed that improved flutter speed prediction 
resulted when measured wing torsional stiffness and sting 
flexibility were incorporated in the mathematical model. 

The excellent nominal and degraded performance of the FSS indi- 
cated that digital implementation of control systems was a viable 
alternative and was capable of supporting multiple advanced 
control concepts. Discrete time and optimal control techniques 
represent the next logical step in the synthesis and implemen- 
tation of active control systems. 
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FIGURE 2-2 - BLOCK DIAGRAM OF FLUHER SUPPRESSION SYSTEM 
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FIGURE 2-4 - REDUNDANT DIGITAL FLUTTER SUPPRESSION SYSTEM SIGNAL PATHS 
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3.0 FLUTTER SUPPRESSION SYSTEM DESIGN CRITERIA AND METHODOLOGY 

This section describes the design criteria and methodology used 
to develop and implement a flutter suppression system capable of 
stabilizing the flutter modes of the B-52E aeroelastic model. 
Paragraph 3.1 presents the criteria and methods used in sel- 
ecting the ballast and sting attachment configuration. Similar 
information is given for the control law synthesis and system 
implementation in Paragraphs 3.2 and 3.3, respectively. 

3.1 Ballast and Attachment Configuration 

The initial effort in the program centered around the structural 
configuration of the model necessary to produce the desired 
flutter characteristics. The constraints imposed on the flutter 
characteristics were as follows: 

• The flutter modes were to exhibit a violent onset, with the 
structural damping slope at flutter to be between approx- 
imately 0.001 and 0.003 per m/s (0.002 and 0.006 per KTAS 
airplane scale). 

t The flutter modes were to have frequencies between 13 Hz 
and 25 Hz. 

• No flutter modes other than the primary wing flutter modes 
(symmetric and antis^metric) were to exist at dynamic pres- 
sures below 4788 N/m^ (100 psf). 

• Wind tunnel testing was to be conducted in 95 percent freon 
with a mass density of 2.58 kg/m^ (0.005 Ib-sec^/ft^) (equiv- 
alent to 6400 m (21 000 feet) atmospheric altitude). 

t In order not to exceed the model design limit dynamic pres- 
sure of 4788 N/m^ (100 psf) the model was to be configured 
to flutter at dynamic pressures below 3831 N/m^ (80 psf). 

Using these constraints the methodology illustrated on Figure 
3-1 was used to define an acceptable configuration. To begin 
the process a start-up configuration was chosen and equations 
of motion were generated. The flutter characteristics of the 
model were then compared to the flutter criteria to determine 
whether the present configuration was acceptable or needed to 
be changed. If a change was indicated the ballast and/or sting 
attachment were updated and the process started over. When the 
flutter criteria were satisfied the structural modifications 
(ballast and/or sting attachment) were designed and transmitted 
to NASA for implementation. Synthesis of the flutter suppres- 
sion system control laws was then ready to begin. 
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Flutter Criteria Met? 


Start-up 

Configuration 



FIGURE 3-1 - CONFIGURATION DESIGN METHODOLOGY 

3.2 Flutter Suppression System Synthesis 

Synthesis of the flutter suppression system (FSS) was centered 
around two tasks, selection of the sensor(s) and control sur- 
face(s) and synthesis of the control laws. The criteria used 
in synthesizing the FSS were as follows: 

• The FSS was to be synthesized in the continuous time domain. 

a The FSS should extend the flutter dynamic pressure at least 
20 percent while maintaining the stability of all other modes 
(both structural and those due to the FSS). 

• The FSS should possess +6 dB of gain margin and ±0.7854 rad 
(45 degrees) of phase margin at and below the flutter velocity. 

0 For digital implementation purposes, the FSS should not have 
modes higher in frequency than approximately 100 Hz. 

With these criteria the methodology illustrated in Figure 3-2 
was used to define the FSS. This methodology consisted of two 
major steps, selection of control surface and sensor parameters 
and synthesis of the control law. The process was started using 
an initial configuration of control surface and sensors which was 
updated until the desired modal coupling was achieved. The next 
step was to synthesize the control laws using an iterative pro- 
cess until the required stability characteristics were obtained. 
During the control law synthesis the modal coupling characteris- 
tics were reviewed with the option of returning to the control 
surface/sensor selection step with revised modal coupling crit- 
eria. After all criteria had been met, the next step was imple- 
mentation of the resulting control law(s). 
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FIGURE 3-2 - FLUHER SUPPRESSION SYSTEM SYNTHESIS METHODOLOGY 

3.3 Flutter Suppression System Implementation 

Implementation of the FSS consisted of designing and implement- 
ing both hardware and software in essentially two parallel, but 
inter-dependent tasks as shown on Figure 3-3. The design crit- 
eria used in this methodology were as follows: 

• The FSS was to be implemented on digital computers using 
difference equations to represent the filters. 

• The frequency response of the FSS should have no significant 
deviation from the ideal around the flutter frequency. 

• The FSS was to have a single fail -operate capability. 

0 The FSS should be capable of detecting when a failure had 
occurred and indicating this fact to the operator. 

0 The sensors and control surfaces were to have no redundancy 
due to size limitations. 

Due to the limitations on the end-to-end redundancy of the 
system the following failures were used in the failure analysis 
as those constituting a single failure. 

0 A single computer failing to update or incorrectly updating 
it's output. 

0 A single failed channel of an A-to-D or D-to-A converter. 

0 A complete failure of an A-to-D or D-to-A converter unit. 

0 A failure in an analog voting device. 

The FSS was required to operate with all but the last failure. 
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FIGURE 3-3 - HARDWARE/SOFTWARE DESIGN METHODOLOGY 



The first step in both the hardware and software design was to 
determine the most accurate and efficient means of implementing 
the FSS. The output of this effort was the information necessary 
to select the analog to discrete filter transform and to define 
the hardware requirements. From this initial step the hardware 
capable of meeting the hardware requirements and the initial 
configuration could be selected. Failure detection analysis was 
then conducted and the hardware configuration modified until all 
failures could be detected. The resulting hardware design was 
then ready to be implemented. 

Design of the software configuration proceeded in parallel with 
the hardware design after the initial requirements were deter- 
mined. After a transform was selected an initial software config- 
uration was chosen and an iterative process was followed until 
the frequency response of the digital filter met the frequency 
response criteria. After a satisfactory digital filter was found, 
the failure detection software was defined using the constraints 
imposed by the hardware configuration. This software, along with 
any support software, was then consolidated and implemented. 
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4.0 AEROELASTIC AND DYNAMIC ANALYSIS 

This section describes the development of structural and 
aerodynamic mathematical models from which equations of 
motion were produced for the B-52E aeroelastic wind tunnel 
model. Flutter analyses were then conducted to define a 
configuration that met the critiera set forth in Section 3.0. 
Structural analysis, including development of the equations 
of motion, is presented in Paragraph 4.1 and flutter analysis 
in Paragraph 4.2. 

4.1 Structural Analysis 

Using the methodology presented in Paragraph 3.1, mathematical 
models were developed that describe the vibrational and aero- 
dynamic characteristics of the wind tunnel model. These math 
models, developed initially in airplane scale, were used to 
produce equations of motion which were reduced to model scale 
using the scale factors in Table 4-1. 

4.1.1 Vibration - Elastic and inertia characteristics of the airplane 
were represented with a lumped parameter idealization. Inertia 
properties (mass, and first and second moments of mass) were 
lumped at the appropriate elastic axis stations. Structural 
stiffness properties were defined by specifying the beam stiff- 
ness parameters El and GJ at each end and the center of each 
beam connecting the elastic axis stations. Tapered beam 
element stiffness representations for the elastic axis were 
generated using the three sets of stiffness properties spec- 
ified for each beam. Cantilevered vibration modes were computed 
for each of the airplane components plus a semi-rigid component 
representing the wind tunnel sting mount. The airplane components 
included the forward and aft body, wing, horizontal stabilizer, 
vertical fin, and wing ballast (when added). The vertical fin 
and horizontal stabilizer were treated as rigid components 

with their mass lumped on the aft body. The forward and aft body 
and wing represented airplane components. 

Coupled vibration modes were determined using a sufficient number 
of canti levered component modes to adequately represent the 
desired low frequency response of the airplane. The equations 
of motion and flutter analyses were based on 32 coupled modes. 

4.1.2 Aerodynamics - Unsteady aerodynamic forces were generated using 
a three-dimensional plate doublet finite element solution. This 
theory accounts for Mach number and finite span effects and includes 
aerodynamic coupling between airplane components. The unknown 
pressure distribution was determined for each airplane mode by 
considering pressure to be a constant over a given aerodynamic 
panel and solving for the pressure based on a specified reduced 
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Scale 


Dimension 


Density 


Froude No. 


Mass Ratio 


Velocity 


Dynamic 

Pressure 


Mach No. 


Frequency 


Weight 


Mass 

Inertia 


Stiffness 


Area 

Inertia 


External 
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Bending 

Moment 


Stress 
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TABLE 4-1 

MODEL SCALE FACTORS 


Formula 


Selected 


Tunnel = .00499 


Airplane Alt, = .0012249 



^ ^ 
Va ^ 


Pm 3 


% /Cm\2 


f^M ^ 


BMm 


Factor 





\ 

^A> 

/ 

pf 

r( 

VMf 

Vf 

_ ®A 


4.07 


1.0 


1.0 


.183 


.136 


.375 


5.48 


151 X 10'® 


.168x 10'® 



1.23 X 10-® 


151 X 10-B 


5.03 X 10’® 


.136 
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frequency and Mach number. The airplane was modeled with 
trapezoidal panels arranged in strips parallel to the free- 
stream. The panel arrangement is shown on Figure 4-1. 

4.1.3 Equations of motion - Initial equations of motion were formed 
using complex oscillatory aerodynamic coefficients generated 
for specific values of the frequency parameter, u/Ug. Final 
equations of motion were formulated in terms of real matrices 
through introduction of an "interpolating" or "approximating" 
function. 


The original equations were the standard form: 

^-(jw)^[Mass] + (joj) [Damping] + [Stiffness] j{q(jti>)} 

+PUq2 [Ai(^)j([Ce]{q(ja.)} +^[Cz]{q.(ju>)} 


where q is the generalized coordinate and Ai is an aerodynamic^ 
influence coefficient matrix which can be evaluated for specific 
values of u/Ug. The matrices, Ce, C^, and Cyj prescribe the 
usual linearized boundary conditions. 

If one of the elements of the complex matrix Aj is plotted, as 
u takes on selected values from 0 to 40 radians/second (airplane 
scale), the result appears as the X's on the sketch below. 


Imaginary Part 


S= 0 + jlO 



S = 0 + j40 
Real Part 


= 0 + j20 


The solid line in the sketch is an approximating function, 
chosen as a rational polynominal function of the complex 
variable S. The circles are values of the approximating 
function at values of S for which the X's are plotted. The 
approximating function was chosen to permit accurate approx- 
imation of the time delays inherent in the unsteady aerodynamics 
subject to the following restrictions: 

t It must have complex conjugate symmetry 

• It must have denominator roots in the left half-plane 
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• It must approximate the value of the complex coefficient when 
S = 0 + jw, for those values of w analyzed. 

The approximating function for each element in the aerodynamic 
influence coefficient matrix was determined after analysis at 
twelve discrete frequencies. When the approximating functions 
are substituted in the equations of motion for the complex 
aerodynamic coefficients, a new set of equations results, whose 
coefficients are coefficients of the approximating function. 

After rearrangement, the final form of the equations of motion 

with variable density p and velocity U and without gust penetration 

is: 


^S^[Mass] + S[Damping] + [Stiffness]j{q(S)} 

4 \ 

+(s2p[Ci] ♦ SPUoCCj] ♦ PU„2[C3] + PU„^ 



The items in the first line of the above equation are structural 
coefficients; items in the second line are aerodynamic coefficients; 
items in the third line are gust velocity coefficients; where: 


S 

P 

Uo 

[Mass] 

[Damping] 

[Stiffness] 

[Ci],[C2],[C3] 

[D] ] »[D2^ »[D3],[D^] 

[Bi],[Gi], 

[Ro].ERi].[R2].[R3].[R43 

q(S) 

Wg(S) 

Vg(S) 


= Laplace variables 

= Air density 

= True airspeed 

= Structural mass 

= Structural damping 

= Structural stiffness 

= Aerodynamic parameters 

= Aerodynamic parameters 

= Lift growth parameters 

= Vertical and lateral gust 
coefficients 

= Rigid body, structural and 
control surface freedoms 

= Vertical Gust 

= Lateral gust 
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Because of the continuity of the aerodynamic coefficients as u 
varies (no aerodynamic poles or zeroes in the vicinity of the 
imaginary axis) these equations are considered to be a good 
approximation of the Laplace transformed equations. They 
should not be depended upon for values of S too remote from 
the imaginary axis or above the highest frequency analyzed 
(100 Hz; model scale). The generalized equations of motion 
were augmented with an additional degree of freedom for each 
control surface. The control surfaces defined for the model 
were inboard, midspan and outboard flaperons, an inboard aileron 
and two outboard ailerons. 

4.2 Flutter Analysis 

Using the equations of motion developed in the previous paragraph, 
analysis was conducted to define a configuration that safisfied 
the flutter criteria. The model modifications (wing ballast and/ 
or sting attachment) were then designed and transmitted to NASA 
for implementation. 

4.2.1 Configuration definition - All analyses were conducted using a 

model scale air density of 2,58 kg/m3 (0,005 Ib-sec^/ft^) 6400 m 
(21 000 feet) equivalent airplane altitude. The initial configur- 
ation had the engine nacelles and external fuel tanks removed from 
the wing and the sting mount located at the existing cable mount 
block. Though this configuration exhibited dual flutter modes 
with violent onset, the flutter speeds occurred well above 4788 
N/m2 (100 psf) design dynamic pressure. Various ballast arrange- 
ments were then Investigated in order to reduce the flutter velo- 
city. Included in the investigation were masses located fore and 
aft of the elastic axis at the outboard nacelle and external fuel 
tank locations. The following ballast arrangement produced satis- 
factory symmetric flutter characteristics. 

t 1.37 kg (3.01578 pounds) ballast attached to the elastic axis 
at the outboard nacelle attach point. 

• 0.05443 kg (0.11974 pounds) ballast attached to the elastic 
axis at the external fuel tank location. 

However, the antisymmetric axis exhibited three flutter modes one 
of which occurred well below the others at 1915 N/m^ (40 psf). 

This 4 Hz mode appeared to be wing chordwise bending coupling 
with fuselage side bending through the moment arm existing between 
the wing attach point (BS 0.05224 (20,5667)) and the sting mount 
(BS 0,7764 (30.5667)). The sting attach point was, therefore, 
moved to the wing attach point which produced satisfactory flutter 
characteristics in both axes. 
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4.2.2 F lutter results - Flutter characteristics were Investigated for 
all dynamic pressures below 4788 N/m^ (100 psf). The frequency 
and damping ratio versus dynamic pressure for both syninetric and 
antisyimietric flutter modes are shown on Figures 4-2 and 4-3. 

Both flutter modes were made up of first wing torsion and second 
wing bending, with the torsion mode becoming unstable at 2873 N/m^ 
(60 psf). No other modes exhibited instability below 4788 N/m'^ 
(100 psf) in either axis. 


4.2.3 Ballast design - The two masses were designed to give the same 

inertia effects as those used in the analysis. Both had circular 
cross-sections to minimize aerodynamic effects and produce maximum 
rigidity. Attachment to the wing was achieved in the same fashion 
as the outboard nacelle and external fuel tank. 
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5.0 FLUTTER SUPPRESSION SYSTEM SYNTHESIS AND ANALYSIS 

This section describes the synthesis and performance analysis 
conducted for the B-52E aeroelastic model FSS. A review of 
the criteria used in the FSS synthesis is given in Paragraph 
5.1. Detailed descriptions of the FSS synthesis and performance 
evaluation are presented in Paragraphs 5.2 and 5.3, respectively. 

5.1 Synthesis Criteria and Constraints 

The FSS was synthesized using the criteria and methodology 
described in Paragraph 3.2. A review of the criteria and 
constraints is presented in the following paragraphs. 

5.1.1 Stability criteria - The stability criteria used in synthesizing 
the FSS were as follows: 

• The FSS should extend the flutter dynamic pressure at least 
20 percent while not significantly degrading the damping of 
any other structural mode. 

• The FSS should possess MIL-F-9490D stability margins below 
the flutter speed; that is +6 dB of gain margin and +0.7854 
rad (45 degrees) of phase margin. 

5.1.2 System constraints - System constraints are usually constraints 
which arise due to physical limitations in the mechanization 
process. The constraints under which the FSS was synthesized 
are given as follows: 

• All modes of the FSS were to be below approximately 100 Hz 
(model scale) to ease digital implementation. 

t The control surface actuation systems were assumed to have 
certain dynamic properties based on experimental results on 
these systems (Reference 5). 

5.2 Flutter Suppression System Synthesis 

This paragraph describes the synthesis of the FSS for the B-52E 
aeroelastic wind tunnel model. Paragraphs 5.2.1 through 5.2.3 
present the control surface and sensor selection and the control 
law synthesis. In Paragraph 5.2.4, the integrated configuration 
of the symmetric and antisymmetric flutter suppression system 
is presented. 

5.2.1 Sensor and control surface selection - Selection of the sensors 
and control surfaces was performed using the methodology des- 
cribed in Paragraph 3.2 and the math models defined in Section 
4.0. Only wing control surfaces, both singularly and in com- 
bination, were investigated, since they are most effective in 
controlling flutter. Besides the existing three segment flaperons 
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and outboard aileron, two additional control surfaces were 
included in both the symmetric and antisymmetric math models. 

As shown on Figure 5-1, the two new control surfaces were an 
aileron located just outboard of the existing one between 
WBL 0.672 (26.47) and WBL 0.716 (28.17) and an inboard 
aileron between WBL 0.343 (13.50) and WBL 0.419 (16.50). 

Also shown on this figure is the sting attach point and the 
two flutter producing ballasts attached at the outboard nacelle 
and external fuel tank locations. 



FIGURE 5-1 - PLAN VIEW OF B-52E AEROELASTIC MODEL 

Since the flutter modes were primarily wing bending and torsion 
with very little fuselage or tail motion, only sensors located 
along the wing were investigated. Though other types of sensors 
were considered, vertical accelerometers were chosen as the 
primary sensor candidates because of size limitations and 
nature of the flutter modes. 

Zero root locus techniques were used to evaluate modal coupling. 
This involved finding the zeroes (numerator root) of the transfer 
function of the sensor response due to a control surface input 
and comparing their locations to those of the open-loop poles 
(denominator roots) on the S-plane. 
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Using this technique, an extensive study was conducted to define 
the most promising combinations of sensors and control surfaces. 

The location of the vertical accelerometers was varied in conjunction 
with the following control surfaces: 

a. The existing outboard ailerons (WBL 0.597 (23.5) to WBL 0.672 
(26.47)) 

b. The new outboard ailerons (WBL 0.672 (26.47) to WBL 0.716 
(28.17)) 

c. The combination of the outboard ailerons (WBL 0.597 (23.5) 
to WBL 0.716 (28.17)) 

d. The outboard segment of the existing three segment flaperons 
(WBL 0.227 (8.92) to WBL 0.297 (11.68)) 

e. The full three segment flaperons (WBL 0.087 (3.42) to WBL 
0.297 (11.68) 

f. The new inboard ailerons (WBL 0.343 (13.5) to WBL 0.415 
(16.33)). 

From this study the most promising combination appeared to be a 
vertical accelerometer at WBL 0.610 (24.0) coupled with the combined 
outboard ailerons (option c above). A zero root locus, illustrating 
this fact is shown on Figure 5-2 where both the combined and existing 
aileron zeroes are plotted as the accelerometer location is varied 
around the nominal position. Although the flutter mode is well 
coupled (by virtue of it's distance from any zeroes), mode qo is 
also well coupled, which is undesirable. By summing the vertical 
acceleration at WBL 0.185 (7.29), with the outboard accelerometer, 
mode qs is shown to be decoupled. 

This study was conducted using the symmetric equations of motion at 
several dynamic pressures with good results in all conditions. The 
antisymmetric equations of motion behaved in essentially the same 
manner allowing the use of the same sensors and control surfaces for 
both systems. Therefore, the selected sensor/control surface con- 
figuration was the sum of the vertical accelerations at WBL 0.610 
(24.0) and WBL 0.185 (7.29) fed back to the combined outboard ailerons. 

The model was modified by NASA to incorporate the larger outboard 
ailerons. The existing outboard ailerons were removed and replaced 
with larger surfaces. From practical structural considerations it 
was necessary to make the new ailerons about 5 percent shorter in the 
spanwise direction than the ailerons used in the analysis. This 
reduction in area, of course, does decrease aileron effectiveness 
slightly, but it should not have a significant effect on the perfor- 
mance of the FSS. The aileron linkages were modified inside the 
fuselage so that each aileron was actuated independently. On the 
original model the ailerons were driven symmetrically by a single 
torque motor. On the modified model the original aileron motor was 
used to actuate the left aileron and what was formerly the right 
flaperon motor was used to actuate the right aileron. 
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5.2.2 Symmetric flutter suppression system synthesis - Synthesis of 
the FSS was conducted using the first (lowest frequency) 18 
degrees of freedom of the math model developed in Section 4,0. 
The open loop flutter characteristics of this model are shown 
on Figure 5.3. Both axes have nearly identical flutter dynamic 
pressures of about 2873 N/m^ (60 psf). Since the goal was to 
extend the flutter dynamic pressure at least 44 percent, synthe- 
sis of the FSS was performed at a point just above flutter, 3112 
N/m^ (65 psf), and just above 1.2 times the flutter dynamic 
pressure, 4309 N/m^ (90 psf). The system was then evaluated 
at the other dynamic pressures as synthesis progressed. 


The synmetric FSS was synthesized using root locis techniques. 
This iterative process involved synthesizing a filter that met 
gain and phase requirements as nearly as possible, evaluating 
it's effects as the feedback gain was varied, and then updating 
the filter equation to improve the results. Included in the 
feedback loop were the dynamics of the aileron actuation system 
as follows: 


^Actual _ 62500 Rad 5_1 

QComraand + 1 50S + 62500 Rad ' 

During the synthesis it was noted that the fairly light damping 
ratio of the actuator mode was causing adverse coupling with modes 
in the same frequency range. This problem was solved by placing 
the following actuator compensation in series with the actuator. 

^ 16096(s 2 + 150S+ (250)^) Rad ^ Eq. 5-2 

^ (S + 250) ($2 + 1700$+ (2006)2) 

This has the effect of making the actuator behave as a first 
order lag at 250 rad/s. 


Using this compensated actuator, the following symmetric F$$ 
filter was derived which satisfied all stability criteria. 


“^Ail Cmd 


W 0.610 (24.0) 

+ WBL 0.185 (7.29) 


(s^) (?^) 

A5Q2 $2 + 12 $ + 4q2\ /$2 + 160$ + 125^ \ 
^ \402 ($ + 150)2 a s2+ 50$ + 1252 / 


where 5 

is positive down 


Cgd 


defined as positive trailing edge down and 
An explanation of the filter terms follow 


5-3 


Rad . 

g 

I 


t The first term in parenthesis is a washout which will remove 
any steady-state commands from the accelerometers. 


• The second and third terms provide 40 dB/ decade of high fre- 
quency gain attenuation. 

t The fourth term provides additional gain and phase lead at 
the flutter frequency. 


• The last term provides a gain peak at the flutter frequency. 
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Figure 5-4 shows the frequency response of the filter and illus- 
trates the effect the last two terms in the filter have on the 
gain at the flutter frequency. Gain root loci of the symmetric 
FSS at 3112 N/m^ (65 psf) and 4309 W (90 psf) illustrating 
the effects of the system on the stability of all modes are 
shown on Figures 5-5 and 5-6. 



Frequency, Hz 



Frequency, Hz 

FIGURE 5-4 - FLUTTER SUPPRESSION SYSTEM FILTER FREQUENCY RESPONSE 
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— I 1 ^ h— • 1- j 300 

-80.0 -60.0 -40.0 -20.0 0 20.0 

a, rad/s 

FIGURE 5-5 - ROOT LOCUS OF THE SYMMETRIC FLUHER SUPPRESSION SYSTEM, 
DYNAMIC PRESSURE = 3112 N/m2 (65 psf) (CONCLUDED) 



Actuator and FSS: Equations 5-1, 5-2 and 5-3 



32 


a, rad/s a, rad/s 

FIGURE 5-6 - ROOT LOCUS OF THE SYMMETRIC FLUTTER SUPPRESSION SYSTEM 
DYNAMIC PRESSURE = 4309 N/m^ (90 psf) 
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5.2.3 Antisymmetric flutter suppression system synthesis - Since the 
basic flutter characteristics of the symmetric and antisymmetric 
axes were virtually identical, it was anticipated that the 
symmetric FSS filter would work for the antisymmetric axis. 

The symneteic system was evaluated on the 18 degree-of -freedom 
antisymmetric math model using the following sign convention: 

• Vertical Acceleration - positive left wing down and right wing 

up. 

• Aileron Displacement - positive left wing trailing edge down 

and right wing trailing edge up. 

The performance was satisfactory as shown by the root loci of 
Figures 5-7 and 5-8. 

5.2.4 System configuration - By performing the appropriate summing 
on the accelerometer signals and actuator commands the FSS can 
be put into the form shown on Figure 5-9. Because the filters 
are identical, a simplified configuration can be obtained which 
treats the flutter modes as left wing and right wing modes 
instead of symmetric and antisymmetric modes, shown on Figure 
5-10. 

5.3 Flutter Suppression System Performance Analysis 

Analysis was conducted to determine the performance of the sym- 
metric and antisymmetric flutter suppression systems. Conducted 
initially at the synthesis conditions of 3112 N/m^ (65 psf) and 
4309 N/m^ (90 psf), this analysis was performed, after the control 
laws were selected, to verify satisfactory performance at all 
conditions below 4788 N/m^ (100 psf). 

5.3.1 Flutter damping performance - The performance of the system in 
damping the flutter mode was evaluated by computing the character- 
istic roots of the model with the FSS operating. This analysis 
also allowed evaluation of the degradation in damping on the 
other structural modes.. 

Plots of damping ratio and frequency of the symmetric and anti- 
symmetric flutter modes with the FSS on and off are given on Figures 
5-11 and 5-12. Note that the flutter speed for both axes has 
been extended beyond 4788 N/m2 (100 psf), with very little 
change in flutter mode frequency. Tables 5-1 and 5-II give 
the damping and freguency of all symmetric and antisymmetric 
modes at 3831 N/m^ (80 psf). Additional data at other dynamic 
pressures is presented in Appendix A. 

5.3.2 System stability margins - During the synthesis process, system 
gain and phase margins were established by use of phase-gain 
root loci. These plots consist of the usual loci associated 
with variations in system gain and additional loci which have 
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been calculated with a given phase shift superimposed over the 
entire S-plane. An example of this type of plot is shown on 
Figure 5-13. 

After the control laws had been selected the gain and phase 
margins for every condition below 4788 N/m^ (100 psf) was estab- 
lished using Bode techniques. This method involves evaluating 
the loop frequency response at the points where the phase is 
3.142 rad (180 degrees) and the gain is 0 dB to determine gain 
and phase margins, respectively. Plots of the minimum gain and 
phase margins for the FSS are shown on Figures 5-14 and 15. 

Note that the margins required at and below 2873 N/m'^ (60 psf) 
are met up to 3352 N/m^ (70 psf). At higher dynamic pressures, 
the FSS was intentionally designed to favor negative phase and 
gain margins. This decision was based on experience with other 
flutter systems where, in general, more phase lag will exist 
than expected and the control surfaces are less effective than 
predicted making it more desirable to have too much gain than 
too little. 

5.3.3 Control surface requirements - Control surface displacement and 
rate requirements were generated using power spectral-density 
(PSD) techniques. Though little is known about the wind tunnel 
turbulence spectrum or amplitude, a rough estimate of the 
control surface requirements was made using a Von Karman spectrum 
with a gust length of 30.48 m (100 feet). By integrating the 
PSD of control surface displacement and rate the RMS values 
were obtained. This data for a 0.3048 m/s (1 ft/sec) RMS 
turbulence level is given on Figures 5-16 and 5-17. Since the 
requirements are fairly constant below 3831 N/m^ (80 psf), an 
indication of the expected control activity above the flutter 
velocity can be obtained at sub-critical speeds. 
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Actuator and FSS: Equations 5-1, 5-2 and 5-3 
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-100.0 -80.0 -60.0 -40.0 -20.0 0 20.0 

a, rad/s 

FIGURE 5-8 - ROOT LOCUS OF THE ANTISYMMETRIC FLUHER SUPPRESSION SYSTEM, 
DYNAMIC PRESSURE = 4309 N/m^ (90 psf) (CONCUDED) 
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Damping Ratio, 5 Frequency, Hz 
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FIGURE 5-n - SYMMETRIC FLUTTER MODE CHARACTERISTICS WITH THE FLUHER SUPPRESSION 
SYSTEM ON AND OFF 
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FIGURE 5-12 - ANTISYMMETRIC FLUTTER MODE CHARACTERIS' 
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TABLE 5-1 

SYMMETRIC MODAL FREQUENCY AND DAMPING 

• Symmetric 

• Dynamic Pressure = 3831 N/m^ (80 psf) 


Mode 

Number 

FSS 

Off 

FSS 

Or 

Frequency, 

Hz 

Damping 

Ratio, 

C 

Frequency, 

Hz 

Damping 

Ratio, 

C 

1 

6.30 

,144 

6,30 

,151 

2 

8.79 

.129 

8,79 

.128 

3 

9.07 

,0025 

9,10 

.0057 

4 

17.51 

.256 

17,64 

.0827 

5 

19.40 

-.0627 

19,57 

.0624 

6 

25.52 

.0169 

25.52 

.0169 

7 

27.50 

.0078 

27.50 

.0078 

8 

28.05 

.0353 

27,17 

.0419 

9 

40.80 

,0093 

40.66 

.0117 

10 

45.20 

.0161 

44.94 

.0194 

11 

45.90 

.0051 

45.90 

.0051 

12 

51.68 

.0373 

53.74 

.0503 

13 

64.54 

.0192 

65.41 

.0164 

14 

72.72 

.0197 

71.78 

.0299 

15 

79.87 

.0069 

79.06 

,0072 

16 

84.94 

.0071 

84.94 

.0071 

17 

99.77 

.0138 

98.55 

.0144 

18 

106.3 

.0076 

106.3 

.0076 

Filter 

— 

— 

36.49 

.255 
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TABLE 5-1 I 

ANTISYMMETRIC MODAL FREQUENCY AND DAMPING 

• Antisymmetric 

• Dynamic Pressure = 3831 N/m^ (80 psf) 


Mode 

Number 

FSS 

Off 

FSS 

On 

Frequency, 

Hz 

Damping 

Ratio, 


Damping 

Ratio, 

? 

1 

6.11 

.154 

6,22 

.172 

2 

7,12 

,105 

7,12 

,105 

3 

9.09 

,0023 

9,11 

,0052 

4 

17,44 

,252 

17,47 

.0693 

5 

17.36 

.109 

17,37 

,111 

6 

19.59 

.0629 

19,79 

.0699 

7 

19.40 

-.0645 

19,59 

.0567 

8 

21.95 

.0078 

21,95 

,0078 

9 

28.07 

.0339 

27,19 

.0443 

10 

32.02 

.0146 

32.02 

.0146 

11 

40.80 

.0092 

40.65 

.0116 

12 

45.21 

.0159 

44.93 

.0191 

13 

51.68 

.0371 

53.47 

.0473 

14 

58.91 

.0114 

58,91 

.0114 

15 

59.99 

.0085 

59,99 

.0085 

16 

64.57 

.0187 

65.32 

.0160 

17 

72.72 

.0191 


.0282 

18 

77.03 

,0075 

77.03 

.0075 

Filter 

— 

— 

35.70 

,250 
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Actuator and FSS: Equations 5-1, 5-2 and:5-3 



FIGURE 5-13 - GAIN/PHASE ROOT LOCUS OF THE SYMMETRIC FLUHER SUPPRESSION SYSTEM 
DYNAMIC PRESSURE = 4309 N/m2 (90 psf) 
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FIGURE 5-14 - MINIMUM GAIN MARGINS OF THE FLUTTER SUPPRESSION SYSTEM 
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Phase Margin, degree 








• Synrietric 

• Gust Length = 30.48m (100 ft.) 
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FIGURE 5-16 - SYMMETRIC CONTROL SURFACE REQUIREMENTS 
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6.0 FLUTTER SUPPRESSION SYSTEM IMPLEMENTATION 

This section describes the mechanization of the FSS control 
laws and fail -operational capability. The control laws, 
which were synthesized in the continuous time domain, were to 
be implemented using discrete time techniques and digital 
computers. The system was to have sufficient redundancy to 
allow non-degraded operation with, and detection of, a single 
failure. 

The overall system requirements and the selected configuration 
are presented in Paragraph 6.1 and the hardware and software 
design are described in Paragraphs 6.2 and 6.3, respectively. 

System performance is presented in Paragraph 6.4 including 
filter frequency response and failure detection. 

6.1 System Configuration 

The methodology presented in Paragraph 3.3 was followed to define 
a configuration which satisfied the design requirements. Initially, 
specific requirements for the hardware were established such as 
computer word length and sample rate. These criteria, together 
with the overall design criteria, were then used to specify the 
hardware and define the FSS configuration. 

6.1.1 System requirements - Some of the more specific configuration 
design criteria, restated from Paragraph 3.3, are as follows: 

• The FSS was to have single fail-operate capability 

• The FSS should be capable of detecting when a failure has 
occurred and indicating this fact to the operator 

• The sensors and control surfaces were to have no redundancy 

Since the entire system could not be made redundant, a single, 
detectable failure was assumed to be one of the following: 

• A single computer failing to update or incorrectly updating 
its output 

• A single failed channel of an A-to-D or D-to-A converter 
t A complete failure of an A-to-D or 0-to-A converter unit 
f A failure in an analog voting device 

A literature survey was conducted to aid in selecting computer 
word length, sample rate, and continuous-to-discrete filter 
transform. As a consequence of this study the following was 
concluded. 
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• The sample rate of the FSS should be at least twice the high- 
est frequency expected at its input (250 Hz) 

0 A computer word length of 16 bits is adequate for systems 
that use scaled integer arithmetic 

0 The best filter transform, in terms of filtering fidelity 
and computation time, is the Bilinear transform (Tustin's 
method). 

A survey of the current methods of implementing redundant control 
systems was also conducted to aid in selecting a redundancy 
scheme that was uncomplicated yet effective. From this study 
came the following requirements: 

0 The system must be triply redundant to produce a fail-operate 
capabi lity 

0 From the reliability and simplicity standpoint, selection of 
one good signal out of three (one possibly bad) is best per- 
formed with an analog device. 

6.1.2 System interface - The FSS interfaces with the wind tunnel model 
via the accelerometer signals and aileron commands. In order to 
provide the proper scaling of these signals and minimize the num- 
ber of A-to-D input channels an analog computer was utilized. 

This interface configuration, shown on Figure 6-1, also provided 
for any pre-or post scaling that was required and an interface 
to aileron excitation generations and data analyzers. 


Aileron Excitation 
!• Data Analyzers 


REDUNDANT DIGITAL — Aileron 
FSS Commands 


• Control Laws 

• Failure Detection/ a 

Indication / 


Analog 
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FIGURE 6-1 - SYSTEM INTERFACE CONFIGURATION 
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Since a considerable number of signals pass between the FSS, ana- 
log computer and analog voter, an interface panel was considered 
a necessity. This panel, shown on Figure 6-2, served as a common 
point at which all the electrical inter-connections were made. 



FIGURE 6-2 - INTERFACE PANEL SIGNAL PATHS 


6.1.3 Final configuration - Several configurations were defined and 

evaluated for their ability to detect the required failures. In 
general, they involved comparing the voter output to the computer 
generated outputs with variations in how the signals were ex- 
changed and cross-compared. The final configuration, which could 
successfully detect all required failures, is shown on Figure 6-3. 

There were two types of signal paths in the final configuration 
’ of the FSS. The FSS filters were implemented by feeding the same 
accelerometer signals through three identical digital filters 
made up of an A-to-D and D-to-A converter and a digital computer. 
The outputs of the filters were then fed through an analog voter 
which produced the aileron commands. Since the voter output would 
always be a good signal unless it failed, it was fed back to each 
of the computers for comparison to other, possibly bad signals. 

The computers then performed the comparison in a circular fashion, 
using the output of the next computer. This "looking-over-each- 
other' s-shoulder" approach allowed detection of any failure with- 
in the FSS itself. A test program would then be used to isolate 
the failed component. 

Using this configuration the hardware required for implementation 
of the FSS was selected. Since the sampling rate was to be at 
least 250 samples per second, the computers had to be capable of 
computing the filtering equations and performing failure detection 
and timing logic in no more than 4 mllli-seconds.' HP 2100 mini- 
computers capable of performing most of its instructions in two 
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FIGURE 6-3 - REDUNDANT DIGITAL FLUTTER SUPPRESSION SYSTEM SIGNAL PATHS 



micro-seconds and an integer multiply in eleven micro-seconds 
were selected. For an average instruction execution time of four 
micro-seconds, this computer could execute 1000 instructions in 
the alloted time, which was more than adequate. Hewlitt-Packard 
A-to-D and D-to-A converters were also used both for interfacing 
ease and because their conversion speeds were higher than required. 

6.2 Hardware Design 

All hardware except the analog voter and the interface panel was 
readily available from laboratory equipment. The design of these 
two components is presented in the following paragraphs. 

6.2.1 Analog voter design - The analog voter was required to select and 
output an unfailed channel of the three input signals. In addition, 
it was to have unity gain at all frequencies and a dynamic range 
of ±10 volts. 

A circuit was selected that produced the middle of the three inputs 
as its output. This circuit, shown on Figure 6-4, works by first 
selecting the maximums of the input signals taken in two's (V/^, 

Vg, Vq). It then outputs the minimum of these three signals which 
would always be the middle of the three inputs. After a success- 
ful breadboard test, three identical channels of the voter circuit, 
along with input/output buffer amplifiers and a switchable inverter 
circuit were mechanized. Circuit diagrams of this circuit and the 
DC power supplies are given in Appendix B. Exterior details of 
the analog voter box are illustrated on Figure 6-5. The front 
panel has main power and on/off switches and input selector knobs 
for each of the three channels. The rear panel contains the FSS 
input, output and test input jacks and provides a mounting place 
for the voltage regulators. The switches for the inverters were 
mounted inside the box to prevent inadvertent contact. 
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6.2.2 Interface panel design - The Interface panel was required to pro- 
vide a common point at which all electrical Interconnections bet- 
ween the redundant computer system, analog voter and analog com- 
puter could be made. The signals routed through the panel were 
symmetric and antisymmetric accelerations, voter inputs (FSS out- 
puts), cross check inputs and voter outputs. The front panel 
details and wiring diagram of the interface panel are illustrated 
on Figure 6-6. The design provided connection to the analog voter 
and analog computer by use of banana jacks on the front panel and 
connection to the digital computer system through terminal strips 
on the rear of the panel. All inter-connection of signals was 
hardwired on the rear of the panel. A bracket to support the 
A-to-D and D-to-A converter cables was also provided. 

6.3 Software Design 

In order to implement the FSS filters on a digital computer the 
describing Laplace domain equations were transformed into differ- 
ence equations. The difference equations along with the failure 
detection and other support software were then programmed on the 
computer in a form that would execute in the reqcfired time. Fin- 
ally, the frequency response and failure detection performance of 
the FSS was assessed. 

6.3.1 Software design requirements - There were two major design re- 
quirements which goverened the overall software configuration. 

The frequency response of the digital filters was required to 
closely match the ideal response in the frequency range of the 
flutter mode. While performing the filter computations the FSS 
was also required to detect any failure within the FSS and give 
some indication to the operator. 

6.3.2 Computer and signal processing equipmeht - The computers chosen to 
Implement the FSS were capable of being -programmed in several high- 
level software languages including FORTRAN and BASIC. However, 

in order to produce the most time-efficient software, a machine 
level language called assembly language was used. This language 
provided direct access to the fundamental instruction set of the 
computer which contained instructions to move data about in memory, 
perform arithmetic and logical operations and enable input-output 
data transfer to peripheral equipment. A brief description of 
the instruction set of the HP 2100 is provided in Appendix B. 

The A-to-D and D-to-A converters were zero-order hold devices, 
that is, the represent the input or output signals in staircase 
fashion. The A-to-D converters were signed, 10-bit devices with 
an input range of ±10 volts and the D-to-A converters were signed, 
12-bit devices with an output range of ±10.24 volts. 

6.3.3 Selection of linear- to- discrete transform - In order to convert 
the Laplace domain equations describing the FSS filter to a 
form that could be programmed on a digital computer a transform 
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WIRING DIAGRAM 

FIGURE 6-6 - INTERFACE PANEL FRONT DETAILS AND WIRING DIAGRAM 




was necessary. The three widely used methods, namely, rectangular 
integration, trapezoidal integration and bilinear transform 
(lustin' s method), were considered in the present study. The 
process of using these three methods to transform a first-order 
lag is illustrated in Table 6-1. In the first two methods the 
Laplace equation was first converted into a block diagram made up 
of summing junctions, gains and integrators. The integration 
approximation was then substituted into the block diagram to give 
the discrete-time approximation ("A is the time delay operator 
where AXp = Xp.-]). From this approximation, the difference 
equation was written. 

In Tustin's method the substitution is the same as for trapezoidal 
integration but it is made directly into the Laplace equation, 
resulting immediately in a difference equation. From this equation 
a block diagram was drawn. 

Tustin's method has certain advantages over the other two tech- 
niques. By comparing the difference equations from Table 6-1 
it can be seen that Tustin's method avoids the delayed feedback 
form which results in the "Yn-2" term in the other equations. This 
"stale data" problem is the primary source of the other method's 
deficiencies. Therefore, Tustin's method was used as the analog- 
to-discrete transform. Difference equations for some common 
filter forms are given in Table 6-II. 

6.3.4 Flutter suppression system filter implementation - Initially, the 
software required to interface with the D-to-A and A-to-D convert- 
ers and to drive a real-time clock was developed. Since all of 
these devices were Hewlitt-Packard equipment, standard interface 
cards and software were available. The real-time clock generates 
interrupts to the computer at a software selectable rate. These 
interrupts were counted in software to determine when the total 
cycle time had elapsed. At the beginning of each cycle the com- 
puter executed the filter equations and failure detection logic 
and counted interrupts until the cycle time had elapsed, where- 
upon the interrupt count was reset and a new cycle begun. 

After successfully transferring analog signals end-to-end through 
the computer system, several single-element filters were program- 
med to gain experience with Tustin's method. Initially, the dif- 
ference equations were programmed using floating point arithmetic 
operations. These filters performed well but consumed a lot of 
time, primarily because floating point arithmetic requires double- 
word length numbers and uses about 20 times more execution time 
than a regular instruction. To overcome this problem the filters 
were programmed using integer arithmetic operations which use 
single-word length numbers and about one-fourth the execution 
time. In order to use this method the input signal and difference 
equation coefficients must be scaled and the software has to check 
for an overflow after some arithmetic operations. 
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TABLE 6-1 

COMPARISON OF DISCRETE TRANSFORMS 




< 

+ 

00 


li 

x|>- 


c 

+ 

oo 


H 

Xl5- 


< 

+ 

OO 


li 


xl>- 


c 

<Ki o 


T3 

o 


3 +J 

c J- 

#d C7^ 
•M CU 
U +j 

<u c 
a: *-H 


« o 

-o •»- 

•r- 4-» 

o 

ISi ^ 

a; CT) 

Q. 0> 

Id 4-> 

u c 


E CO 

X- S- - ^ 
O C T3 
0> O 

C CO 4-) ^ 
c t/i 4J 
»— <0 3 OJ 
S- I- Z 
CO I— ^ 


62 



TABLE 6- I I 

DIFFERENCE EQUATIONS FOR COMMON FILTER ELEMENTS 


Filter 

Difference Equation 

LAG 

A/(S + A) 

Yn = KiYn_i + K 2 (X, + X,.i) 
Ki = (2 -AT)/(2 + AT) 
K 2 = AT/(2 + AT) 

WASHOUT 
S/(S + A) 

Yn = KiY„.i + K 2 (X„ + X^_i) 
Ki = (2 - AT)/(2 + AT) 
K 2 = 2/(2 + AT) 

LEAD-LAG 
(S + B)/(S + A) 

Yn = KiYn_l + K 2 Xn + K3Xn_i 
Ki = (2 - AT)/(2 + AT) 
K 2 = (2 + BT)/(2 + AT) 
K 3 = (-2 + BT)/(2 + AT) 

SECOND ORDER 
a 2 

S2 + 2ABS + a2 

Yn = KlYn_i + K 2 Yn -2 + K3(Xr, + 2 Xn.i + X,. 2 ) , 

Ki = ( 8 - 2 a 2 t 2 )/D 
K 2 = (-4 + 4ABT-a2t2)/D 
K 3 = A^T^/D 
D = 4 + 4ABT + A^t2 

NOTCH 

S^ + 2ACS + a2 
S2 + 2ABS + a2 

Yn = KlY„.i + K 2 Y ,.2 + K3X, + K4Xn_i + K5X,.2 
Ki = ( 8 - 2 a 2 t 2 )/D 
K 2 = (-4 + 4ABT-a2t2)/D 
K 3 = {4 + 4ACT + a2t^)/D 
K 4 = {- 8 + 2 a 2 t 2 )/D 
K 5 = (4-4ACT + a 2 t 2 )/D 
D = 4 + 4ABT + a2t2 















The FSS filters were programmed using the scaled integer techni- 
que and the difference equations found in Table 6-III. All the 
filter coefficients with magnitude less than one were scaled the 
maximum amount, 2l5. Those coefficients greater in magnitude were 
scaled by 214 and the factor of two added after the multiplication 
by shifting the number left one bit. The result of each multipli- 
cation was a 32-bit word which was shifted left one bit to correct 
the scaling and then truncated to a 16-bit word. The filters were 
implemented in a cascade or series fashion with the output of one 
filter serving as the input to the next. 

After the filters were programmed the end-to-end frequency response 
was evaluated. Initially, considerable effort was expended in ad- 
justing the intermediate gains of the filter to eliminate overflow 
conditions. Even after all overflows had been eliminated the fre- 
quency response of the filter deviated significantly from the ideal, 
especially near the flutter frequency. The problem was found to be 
the limited range of the input device (10-bit A-to-D) and the wide 
difference in the frequency responses of the filter terms. The 
first-order lag at 10 rad/s has a very low gain at the flutter 
frequency of 125 radians/second. This reduces the resolution of 
the inputs from 10 bits to about six. When this signal was passed 
through the inverse notch at 125 rad/s, the round-off error dis- 
torted the output considerably. Changing the order of computation 
of the filter terms did not improve the result. 

Since the series implementation would not work an alternate con- 
figuration was devised. This consisted of a parallel implement- 
ation of the filter by performing a partial fraction expansion on 
the filter. The resulting configuration is shown on Figure 6-7. 

The gains prior to the final summation were scaled by 70 and the 
O-to-A gain to reduce their magnitudes to less than one. The 
gain of 70 could then be moved to the analog computer and the 
internal gains scaled by 215 for integer multiplication, as before. 
The multiple terms in the denominator of the original filter 
produced a partial fraction expansion with only three unique 
terms. Since in Tustin's method overflow detection is riot 
required in first-order lags, the only overflow detection was 
performed internal to the second-order term and at the output 
of the final sunmation. The difference equations to implement 
this filter are given in Table 6-IV. 

These filter equations were programmed using a sample rate of 
500/s and the system frequency response evaluated. The phase 
and gain matched the ideal responses closely except at the higher 
frequencies. The execution time delay caused the phase to lag 
the desired phase and the sampling rate caused the gain to go to 
zero at the nyquist frequency (250 Hz). A non-recursive approx- 
imation to a time advance was added in series with the output of 
the filter as follows: 

Yn = 1.875Xn - 1.25Xn_i + 0.375Xn_2. 
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TABLE 6-1 n 

DIFFERENCE EQUATIONS FOR THE FLUHER SUPPRESSION SYSTEM FILTER 


Filter 


Difference Equation (T=.004) 


S 





10 





150 





22500 $2 + 12S + 1600 


1600 $2 + 300S + 22500 


$2 + 160$ + 15625 


S2 + 50S + 15625 


= KlYn-l + K2(Xn-Xn-l) 

Ki = 0.960784(2^^) = 31483 

K2 = 0.980392(2-5) = 32125 


= KiYn.i + K2(Xn + Xn.l) 

Ki = 0.960784(215) = 31483 
K 2 = 0.019078(215) = 643 


= KiYn_i + K2(X„+Xn_i) 

Ki = 0.538462(215) = 17644 
K 2 = 0.230769(215) = 7562 


= G [KiYn.i + K 2 Yn .2 + K3Xn + K4Xn-l + K5Xn.]] 
G = 14.0625 

Ki = 1.0769(21^) = 17643 
K 2 = -0.289941(215) = -9468 
K 3 = 0.609704(215) = 19978 
K 4 = -1.175858(21^) = -19625 
K 5 = 0.58130(215) = 19647 


' KiY„_ 1 + K2Yp_2 + K3Xn + K4Xn_i + KgX^.z 
Ki = 1.6129(21^) = 26426 
K2 = -0.827957(215) = -27130 
K3 = 1.18925(21^) = 19485 
K4 = -1.6129(2^^) = -26426 
K5 = 0.63879(215) = 20929 
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FIGURE 6-7 - PARALLEL EXPANSION- OF THE FLUTTER SUPPRESSION SYSTEM FILTER 


TABLE 6 -IV 

DIFFERENCE EQUATIONS FOR THE PARALLEL FLUTTER SUPPRESSION SYSTEM FILTER 


Filter 

Difference Equations (T= 0,002) 

10 

S+ 10 

Yn = KiYn-l + K2(Xn + Xn-l) 

Ki = 0.980198(215) = 32 H 8 
K 2 = 0.009901(215) = 325 

150 

S + 150 

Yn = KlYn-l + K 2 (Xn + Xn.i) 

Kj = 0.73913(2^^) = 24220 
K 2 = 0.130435(21^) = 4274 

34.3825(5 + 56.811) 
(S2 + 50S + 15625) 

Yn “ KiY^_l + K 2 Yn _2 + K3Xn + K4Xn_i + K5Xn_2 
Ki = 1.84751(21^) , 30270 
K 2 = -0.90616(21^) = .29693 
K 3 = 0,27276(215) = 8938 
K 4 = 0.029323(215) = 951 
Kg s -0.243434(21^) = -7977 
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This prediction algorithm caused the output to be approximately 
what it would have been had no time elapsed between input and 
output. This improved the frequency response to an acceptable 
level . 

In the parallel expansion of a filter a change in any term of the 
filter will affect the entire parallel filter. This prompted the 
development of an alternate parallel filter which did not include 
the washout or the first-order lag at 150 rad/s. These terms 
would then be mechanized on the analog computer for pre- or post 
filtering should this be necessary during the wind tunnel test. 
The configuration of this parallel filter is shown on Figure 
6-8 and the difference equations appear in Table 6-V. 


ACCELERATION 

INPUT 



Aileron 

'Conmand 


FIGURE 6-8 - PARALLEL EXPANSION OF THE PARTIAL FLUTTER SUPPRESSION SYSTEM 
FILTER 


6.3.5 Failure detection and indication - The fail-operate and failure 
detection capabilities of the FSS were primarily determined by 
the hardware configuration as described in Paragraph 6.1. 

However, the detection and indication of a failure was performed 
in software. This software was then integrated with the filter 
software to complete the software design. 

The failure detection software was programmed to detect two 
basic types of failures, erroneous system output and execution 
time overruns. The first type of failure was detected by com- 
paring the absolute difference between the voter and computer 
outputs to a pre-determined threshold. If the error exceeded the 
threshold a failure was declared. Execution time overruns 
were detected by requiring the computer to set a flag at the 
end of execution of the program. If, at the start of the next 
execution interval this flag was not set, a failure was declared. 

The occurrence of a failure or failures was indicated by flashing 
the computer's front panel lights. In order to make the flashing 
rate independent of the failure rate, failures were accumulated 
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TABLE 6 -V 

DIFFERENCE EQUATIONS FOR THE PARTIAL PARALLEL FLUTTER SUPPRESSION SYSTEM 

FILTER 


Filter 

Difference Equation (T = 0,002) 

10 

Vn = KiY„_i + K 2 (X„tX„.l) 

S + 10 

Ki = 0,980198(2^5) = 32118 
K 2 = 0.009901(215) = 325 

150 

Yn = KiYn_i + K 2 (Xn + Xn-l) 

S + 150 

Ki = 0.73913(215) = 24220 
K 2 = 0.130435(215) = 4274 

43.68(5 - 32.938) 

Yp ” l^lYp-l *^ 2 Yn -2 *^3^n ^4^n-l *^5^n-2 

(S2 + 50S + 15625) 

Ki = 1.8475(21^) = 30270 
K 2 = -0.90616(215) = -29693 
K 3 = 0.31712(215) = 10391 
K 4 = -0.02160(215) = -708 
K 5 = -0.33872(21^) ^ .11099 


for about one-half second. Unless no failures had occurred dur- 
ing that time interval the front panel lights were flashed at 
1 Hz. This failure summing technique prevented the flash rate 
from being too fast or slow to be seen by the operator. 

A separate program was written which allowed easy testing of 
each computer system to determine the origin of a failure. List- 
ings of this program and the main program are given in Appendix 

B. 

6.4 System Performance 

The performance of the entire system was evaluated to verify that 
all performance criteria had been met. The testing included the 
frequency response of the system with and without failures, 
analog voter performance and the successful detection of all 
required failures. 
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6.4.1 Filter frequency response - The frequency response of the FSS 
filters was evaluated using a digital transfer function analyzer. 
The resulting plots matched the theoretical response closely as 
shown on Figures 6-9 and 6-10. The deviations at the higher 
frequencies are attributable to effects of the sampling rate on 
phase and the output prediction algorithm on gain. 

6.4.2 Analog voter performance - Initially analog voter performance 
was tested using a frequency response analyzer. Frequency 
responses of the analog voter with and without a failed chan- 
nel showed unity gain characteristics up to 1000 Hz. In order 
to test the operation of the system with one channel failed, 
end-to-end frequency responses of the FSS filters were run. 

The resulting plots, shown on Figures 6-11 and 6-12, showed very 
little degradation due to a single failure. The variations in 
gain at 6.5 Hz and 100 Hz are attributable to the low amplitude 
of the output coupled with variations in voter electronics. 

6.4.3 Failure detection performance - The ability of the FSS to detect 
internal failures was tested by introducing failures within the 
system and observing the results. The following failures were 
introducted by breaking wiring connections and switching off 
components : 

• Failure of an A-to-D or D-to-A converter channel or the 
entire unit 

• Failure of a digital computer 

• Failure of the analog voter. 

In all cases, with an input of reasonable magnitude, the FSS 
was able to detect and indicate to the operator that a failure 
had occurred. 
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Input Scaling - 2.5 V/g 

Output Scaling - 0.02618 rad/V (1.5 deg/V) 



Frequency, Hz 

FIGURE 6-9 - FREQUENCY RESPONSE OF SYMMETRIC FLUHER SUPPRESSION SYSTEM 










Input Scaling - 2.5 V/g 

Output Scaling - 0.02618 rad/V (1.5 deg/V) 



Frequency, Hz 

FIGURE 6-10 - FREQUENCY RESPONSE OF ANTISYMMETRIC FLUTTER SUPPRESSION SYSTEM 
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7.0 


TEST SUPPORT AND RESULTS 


This section describes the wind tunnel test of the FSS at NASA 
Langley Research Center. Support was provided for the initial 
setup and checkout and operation of the FSS during the wind 
tunnel test. Paragraph 7.1 describes the initial setup of the 
FSS prior to the test and the modifications made during pre- 
liminary testing in the wind tunnel. Results of the test are 
presented in Paragraph 7.2 and of the post-test analysis in 
Paragraph 7.3. 

7.1 Flutter Suppression System Preparation and Model Modifications 

Prior to the actual wind tunnel test the wind tunnel model and 
FSS hardware were integrated and tested. This allowed sub- 
system interfaces to be established, wiring to be completed 
and an end-to-end check to be performed. Included in the 
pretest activities was the programming of the analog computer 
and setup of the aileron actuation system. During preliminary 
testing several modifications were also identified which would 
provide a better correlation between analytical and actual 
model characteristics. 

7.1.1 Flutter suppression system preparation - Upon arrival at the test 
site, the FSS was complete with the exception of programming the 
analog computer and readying the aileron actuation system. The 
following functions were programmed on an EAI 580 analog computer: 

• Aileron actuator compensation 

t Accelerometer scaling and summing to form symmetric and anti- 
symmetric signals 

• FSS pre- or post filtering 

• Aileron command scaling and summing to form left and right 
commands 

t Interface to aileron excitation generators (sweeps, steps) and 
data analyzers. 

A patching diagram of the final configuration is provided in Ap- 
pendix B. 

The aileron actuation system is an electro-mechanical system made 
up of electric servo motors and position and rate transducers 
mechanically linked to the control surface (Reference 2). Posi- 
tion and rate feedback loops are closed through electronics loc- 
ated in the tunnel control room. The feedback gains were adjucted 
to give the desired command sensitivity and dynamics. The actuator 
compensation on the analog computer was then adjusted to cancel 
the actuator mode. 
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During initial integration of the model and the FSS it was found 
that the accelerometer outputs contained high frequency noise and 
had a strong DC drift with temperature. Since the input to a dig- 
ital filter needs to be band-limited and centered in the input 
range of the A-to-D converter the alternate FSS was used during 
the wind tunel test. 

7,1.2 Model Modification - During preliminary wind tunnel testing two 

modifications were made to the model to produce closer correlation 
with the associated math models used during analysis. During 
initial setup the right aileron actuation system displayed 
considerably more inertia than the left. This caused large 
changes in actuator dynamics when the torque limit of the servo- 
motor was reached. During preliminary testing the wind tunnel 
turbulence produced large enough aileron commands to cause an 
instability. The aileron system was modified to reduce the 
inertia by re-routing and shortening the mechanical linkage 
between the servo-motor and the control surface. This resulted 
in an inertia roughly equivalent to the left aileron system. 

After this modification both acuation systems were tested and 
found to be capable of full deflections at the flutter frequency. 
The final actuator compensation was as follows: 

. . _ 14348(5^ + 200$ + 95000) Rad 

^ ' (S + 350)2(s2 + 17005 + 4.0 X 106) Rad 

As testing resumed it became apparent that flutter would occur 
at a higher dynamic pressure and model frequency than had been 
predicted by structural analysis. A comparison of GVT and 
analytical data (Table 7-1) showed that the analytical flutter 
pair was lower in frequency than the actual model. Because 
the mass properties of the model were well known, the torsional 
stiffness of the model was assumed to be the source of the 
problem. Since this could not be modified directly,' the frequency 
was lowered to the analytical value by lowering the inboard 
flutter ballast 0.01524m (0.6 inches). This modification pro- 
duced structural frequencies very close to analytical values 
but did not lower the flutter velocity. 

TABLE 7-1 

COMPARISON OF GVT AND ANALYTICAL FREQUENCIES 


Symmetric Mode 

Frequency, Hz 

GVT 

Analysis 

First Wing Bending 

3.3 

3.3 

Second Wing Bending* 

13.5 

12.7 

First Wing Torsion* 

25.2 

23.3 


♦Flutter Pair 
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7.2 


Test Results 


The primary goal of the wind tunnel test was to demonstrate the 
capability of the FSS to suppress dual flutter modes (synmetric 
and antisyiranetric) with violent, nearly-simultaneous onset, and, 
furthermore, that this system could be successfully implemented 
using digital computers as the feedback filters. A secondary 
goal was to investigate the characteristics of the redundant 
system under degraded system operation. 

7.2.1 Flutter mode damping performance - The performance of the FSS in 
damping the flutter modes was evaluated initially by observing 
its effect just below flutter speed. This was done primarily to 
establish the open loop flutter speed and to gain confidence that 
the system was capable of stabilizing the flutter modes. Although 
the flutter dynamic pressure had been inaccurately predicted, 2873 
N/m2 (60 psf) instead of 3926 N/m^ (82 psf), the FSS demonstrated 
good damping capability as shown on Figure 7-1. This strip-chart 
recording shows the response of the flutter mode building up with 
the FSS off and then becoming highly stable when it is turned on. 
This fact is also demonstrated on Figure 7-2 where the symmetric 
acceleration frequency response from an aileron sweep is shown. 
With the FSS off the response of the flutter mode at 19 Hz domin- 
ates the plot but is virtually eliminated when the system is turned 
on. 

After verifying that the FSS was operating properly the dynamic 
pressure of the wind tunnel was increased to values above the 
flutter velocity. As speed was increased the flutter mode remain- 
ed stable with some increase in control surface activity due to 
tunnel turbulence. At the structural limit of the model (4788 N/m^ 
(100 psf)) the control surface activity was approximately 0.0873 
rad (5 degrees) RMS with 0.1745 rad (10 degrees) peaks. Also, a 
mode at about 48 Hz was becoming increasingly active as the tunnel 
speed was increased. It appeared that this mode, probably either 
a structural mode excited by the actuator mode or the actuator 
mode itself, would have defined the system-on flutter boundary had 
the model's structural limit been higher. The maximum dynamic 
pressure tested was 4884 N/m^ (102 psf) which represented an in- 
crease of 24 percent over the FSS off flutter dynamic pressure 
of 3926 N/m2 (82 psf). 

7.2.2 Degraded system performance - The performance of the system was 
evaluated while various degradations were introduced into the 
redundant portion of the FSS. Induced degradations included a 
single channel failure and gain reductions and phase changes in 
one channel while another was failed to maximum input level. The 
effects of these degradations on the output of the voter with sin- 
usoidal inputs are illustrated in Table 7-II. The phase shift 
degradation was introduced by selecting the break frequency of a 
first-order lag in series with the FSS to give the desired phase 
change with unity gain at the flutter frequency (~20 Hz). The 
tests were conducted at 3831 N/m^ (80 psf) where the model was 
stable but the effects of the FSS were obvious. 
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FIGURE 7-1 - MODEL RESPONSE WITH FLUHER SUPPRESSION SYSTEM ON AND OFF 
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TABLE 7-II 

VOTER OUTPUTS WITH DEGRADED INPUTS 
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Failing a single channel had no effect on the performance of the 
FSS as expected. The failed channel was tied to a high level 
(10 volts) and the gain on a second channel incrementally reduced. 
This had the effect of reducing the total FSS gain by half of 
the value of the single channel. That channel was then restored 
to full gain and a series of first-order lags with the appropriate 
gain inserted. This appeared to have the effect of changing the 
phase of the total FSS by half of this amount. These results 
indicated that first the FSS has fail-operate capability and that 
second, with one channel failed high, the FSS would be fairly 
insensitive to degradation in another channel. 

7.3 Post Test Analysis 

Analysis was conducted after the wind tunnel test to determine 
why the flutter speed was higher than predicted. The analytical 
performance of the improved math model was then compared to wind 
tunnel results. 

7.3.1 Changes to structural model - As was noted earlier the difference 
in structural frequencies of the model and math model were 
probably due to a difference in torsional stiffness. Also 

the GVT data indicated that there was some flexibility in the 
sting mount, which had not been modeled. Two math models were 
developed to investigate the effects of, and sensitivity to, 
adding torsional stiffness and/or sting flexibility. The first 
model had only torsional stiffness added to give the same first 
torsion mode frequency. A second math model was developed that 
had some sting flexibility added and the torsional stiffness 
adjusted to math model frequencies. 

7.3.2 Comparison with test results - The math model with torsional 
stiffness added had a flutter dynamic pressure of 3735 N/m^ 

(78 psf) only 3 percent low in velocity. The flutter character- 
istics were similar to the model as illustrated on Figure 7-3. 

This demonstrated the sensitivity of the model to torsional 
stiffness . 

The second math model fluttered at about 4788 N/m^ (100 psf), well 
above the test results. Since the sting flexibility could only 
be guessed, the results were used as an indication of sensitivity 
to this parameter. 

Using the math model with increased torsional stiffness, closed 
loop analysis was conducted to verify the performance of the FSS. 
As shown by the root locus on Figure 7-4, the flutter mode is 
driven very close to its zero when the FSS is turned on. This 
would tend to reduce its acceleration response to an aileron 
sweep dramatically, which was precisely what occurred in the 
actual test. The performance of the FSS in damping the flutter 
mode is illustrated on Figure 7-5. The system-on flutter speed 
is greater than 4788 N/m'^ (100 psf) as was the case in the wind 
tunnel test. 
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FIGURE 7-3 - SYMMETRIC FLUTTER CHARACTERISTICS - MODIFIED STIFFNESS 



Actuator and FSS: Equations 5-1, 5-2 and 5-3 



CT, rad/s a, rad/s 

FIGURE 7-4 - ROOT LOCUS OF FLUHER SUPPRESSION SYSTEM USING UPDATED EQUATIONS, 
DYNAMIC PRESSURE = 3831 N/m^ (80 psf) 
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FIGURE 7-4 - ROOT LOCUS OF FLUTTER SUPPRESSION SYSTEM USING UPDATED EQUATIONS, 
DYNAMIC PRESSURE = 3831 N/m2 (80 psf) (Concluded) 
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FIGURE 7-5 - FLUTTER SUPPRESSION SYSTEM SYSTEM FLUTTER MODE DAMPING 
PERFORMANCE 


In short, the structural analysis indicated that the math model 
was sensitive to both torsional stiffness and sting flexibility. 
The closed loop analysis showed several points of correlation 
to the actual test results. 


8.0 CONCLUSIONS AND RECOMMENDATIONS 


8.1 Conclusions 

Major conclusions resulting from this study are listed below: 

1. The FSS was successful in stabilizing symmetric and antisym- 
metric flutter modes which exhibited violent onset and identical 
flutter velocities even though the experimental flutter velocity 
was considerably different than predicted analytically. This can 
be attributed to the wide stability margins of the FSS. Post 
wind tunnel test analysis indicated that the flutter speed of 
the mathematical model was sensitive to wing torsional stiffness 
and sting flexibility. Improved flutter velocity predictions 
were obtained when measured wing torsional stiffness and sting 
flexibility were incorporated in the model. 

2. Digital implementation of control system filters which were 
synthesized using classical techniques can be performed with good 
frequency response fidelity, especially when a prediction algo- 
rithm is used to reduce time delay effects. 

3. A triple redundant configuration which utilizes a circular 
failure detection scheme produces a fail-operate system capable 
of detecting any single internal failure with indication to the 
operator. 

4. Using the parallel expansion technique in implementing digital 
filters eliminates the need to adjust internal gains to prevent 
under or overflows and thus reduces the total time needed to 
implement the filter. 

8.2 Recommendations 

The recomnendations listed below are offered to suggest areas of 
future research and to ensure the success of these projects. 

1, When testing flutter suppression systems every attempt should 
be made to attain an accurate mathematical model, particularly in 
regard to wing stiffness and sting flexibility. 

2, Stability margins on future control systems should be as wide 
as possible in order to produce systems that are insensitive to 
variations in mathematical models. 

3, The performance and flexibility of this system indicate that 
other ACT (Active Controls Technology) concepts such as gust and 
maneuver load alleviation and relaxed static stability could be 
synthesized using classical control techniques and implemented 
using digital computers. 
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4. As the technology advances, discrete time and advance control 
concepts (Z-transform, optimal control) should be used in the 
synthesis of control systems as well as in their implementation. 
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APPENDIX A 

FLUTTER SUPPRESSION. SYSTEM PERFORMANCE 

This appendix contains additional performance data for the FSS using the 
original equations of motion. Modal damping and frequencies with the FSS 
on and off appear in Tables A-I to A-VI. 

TABLE A-I 

MODAL DAMPING AND FREQUENCIES 


• Symmetric 

• Dynamic Pressure = 0 N/m^ (0 psf) 


Mode 

Number 

FSS Off 

FSS 

On 

Frequency. 

Hz 

Damping 

Ratio. 

Frequency, 

Hz 

Damping 

Ratio, 

C 

1 

3.34 

.0050 

3.34 

.0049 

2 

6.94 

.0075 

6.94 

,0075 

3 

9.22 

.0050 

9.22 

.0050 

4 

12.70 

.0049 

12.70 

.0049 

5 

23.35 

.0050 

23.35 

.0051 

6 

25.38 

.0075 

25.38 

.0075 

7 

27.69 

.0075 

27.69 

.0075 

8 

28.34 

.0050 

28.34 

.0050 

9 

41.14 

.0050 

41.14 

.0050 

10 

45.84 

.0050 

45.84 

.0050 

11 

46.74 

.0075 

46.74 

.0075 

12 

50.93 

.0050 

50.89 

.0049 

13 

64.96 

.0050 

64.92 

.0051 

14 

72.85 

.0050 

72.85 

.0050 

15 

79.85 

,0050 

79.89 

.0048 

16 

85.14 

.0075 

85.14 

.0075 

17 

100.05 

.0050 

100.15 

.0042 

18 

106,3 

.0075 

106.3 

.0075 

Filter 

— 

— 

19.90 

.2000 
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TABLE A-II 

MODAL DAMPING AND FREQUENCIES 


• Symmetric 

t Dynamic Pressure = 1436 N/m^ (30 psf) 


Mode 

Number 

FSS 

Off 

FSS On 

Frequency, 

Hz 

Damping 

Ratio, 

C 

Frequency, 

Hz 

Dampi ng 
Ratio, 
C 

1 

4.56 

.126 

4.49 

,285 

2 

7.70 

.0904 

7.70 

.0904 

3 

9.18 

.0053 

9.19 

.0055 

4 

14.18 

.0935 

14.83 

,0907 

5 

22.18 

.0200 

22.18 

,0194 

6 

25.43 

.0124 

25.43 

.0124 

7 

27.62 

.0077 

27.62 


8 

28.25 

.0261 

28.10 

,0323 

9 

41.04 

.0075 

41.02 

,0076 

10 

45.64 

.0129 

45.56 

.0129 

11 

46.48 

.0063 

46.48 

.0063 

12 

51.22 

.0238 

51,74 

.0288 

13 

64.82 

.0131 

65.04 

,0l39 

14 

72.83 

.0134 

72.62 

.0152 

15 

79.86 

.0061 

79.59 

.0054 

16 

85.06 

.0075 

85.06 

.0075 

17 

99.97 

.0103 

99.51 

.0089 

18 

106.3 

.0075 

106.3 

.0075 

Filter 

— 

— 

23.39 

.338 
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TABLE A-III 

MODAL DAMPING AND FREQUENCIES 


• Symmetric 

• Dynamic Pressure = 2873 N/m^ (60 psf) 


Mode 

Number 

FSS 

Off 

FSS On 

Frequency, 

Hz 

Damping 

Ratio, 

Frequency, 

Hz 


1 

5.61 

.140 

5,82 

.207 

2 

8.37 

.117 

8,36 

.117 

3 

9.13 

,0044 

9.14 

.0058 

4 

16.29 

.163 

16,55 

.103 

5 

20.27 

.0015 

20.74 

,0301 

6 

25.48 

.0151 

25.48 

,0151 

7 

27.55 

.0077 

27,55 

.0077 

8 

28,12 

.0329 

27,51 

,0435 

9 

40.90 

.0087 

40,83 

,0096 

10 

45.38 

.0153 

45.20 

.0168 

11 

46.14 

.0057 

46,14 

.0058 

12 

51 . 49 

.0326 

52,82 

.0422 

13 

64.65 

,0171 

65,21 

.0164 

14 

72.77 

.0175 

72.19 

.0233 

15 

79.87 

.0066 

79.28 

.0065 

16 

84.98 

.0073 

84.98 

.0073 

17 

99.85 

.0126 

98.93 

.0123 

18 

106.3 

.0076 

106.3 

.0076 

Filter 

— 

— 

32.91 

.326 
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TABLE A-IV 

MODAL DAMPING AND FREQUENCIES 


• Antisymnetric 

f Dynamic Pressure = 0 N/m^ (0 psf) 


Mode 

Number 

FSS 

Off 

FSS On 

Frequency, 

Hz 

Damping 

Ratio, 

? 

Frequency, 

Hz 

Damping 

Ratio, 

C 

1 

3.34 

.0050 

3.34 

.0049 

2 

5.30 

.0075 

5.30 

.0075 

3 

9.22 

.0050 

9.22 

.0050 

4 

12.70 

.0049 

12.70 

.0059 

5 

15.94 

.0075 

15.94 

.0074 

6 

21.21 

.0075 

21.21 

.0074 

7 

22.39 

.0075 

22.39 

.0075 

8 

23.34 

.0050 

23.34 

.0051 

9 

28.34 

.0050 

28.34 

.0050 

10 

32.17 

.0075 

32.17 

.0075 

11 

41.13 

.0050 

41.13 

.0050 

12 

45,84 

.0050 

45.84 

.0050 

13 

50.93 

.0050 

50.89 

.0049 

14 

58.84 

.0075 

58.84 

.0075 

15 

59,97 

,0075 

59.97 

.0075 

16 

64.95 

.0050 

64.92 

.0051 

17 

72.86 

.0050 

72.85 

.0050 

18 

77.11 

.0075 

77.11 

.0075 

Filter 

— 

— 

19,88 

.198 

1 
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TABLE A-V 

MODAL DAMPING AND FREQUENCIES 
• Antisymmetric 

§ Dynamic Pressure = 1436 N/m^ (30 psf) 



FSS 

Off 

FSS On 

Frequency, 

Hz 

Damping 

Ratio, 

C 

Frequency, 

Hz 

Damping 

Ratio, 

C 

1 

4.49 

.132 

4.38 

.298 

2 

6.05 

.0827 

6.05 

.0827 

3 

9.18 

.0051 

9.19 

.0054 

4 

14.17 

.0926 

14,84 

.0885 

5 


.0583 

16.33 

.0583 

6 


.0495 

20.77 

.0495 

7 

22.17 

.0194 

22.18 

.0188 

8 

22.23 

.0075 

22.23 

.0075 

9 

28.25 

.0257 

28.12 

.0321 

10 

32.13 

.0116 

32.13 

.0116 

11 

41.04 

.0074 

41.01 

.0076 

12 

45.64 

.0127 

45.56 

.0128 

13 

51.22 

.0237 

51.62 

.0284 

14 

58.88 

.0098 

58.88 

.0098 

15 

59.98 

.0081 

59.98 

.0081 

16 

64.83 

.0129 

65.05 

.0135 

17 

72.84 

.0132 

72.65 

.0152 

18 

77.10 

.0074 

77.10 

.0074 

Filter 

— 

— 

23.39 

.334 















TABLE A-VI 

MODAL DAMPING AND FREQUENCIES 


• Antisymmetric 

t Dynamic Pressure = 2873 N/m^ (60 psf) 


Mode 

Number 

FSS 

Off 

FSS 

On 

Frequency, 

Hz 

Damping 

Ratio, 

C 

Frequency, 

Hz 

Damping 

Ratio, 

C 

I 

5.48 

.149 

5.76 

.230 

2 

6.72 

.0995 

6.72 

.0995 

3 

9.13 

.0041 

9.14 

.0055 

4 

16.25 

.161 

16.54 

.0963 

5 

16.85 

.0860 

16.85 

.0862 

6 

20.15 

.0629 

20.15 

.0627 

7 

20.26 

-.0004 

20.78 

.0296 

8 

22.07 


22.07 

.0076 

9 

28.13 

.0319 

27.56 

.0449 

10 

32.06 

.0135 

32.06 

.0135 

11 

40,90 


40,82 

,0096 

12 

45.38 

.0151 

45.19 

.0167 

13 

51.49 

.0325 

52.70 

.0405 

14 

58.90 


58.90 

.0109 

15 

59.98 

.0084 

59.98 

.0084 

, 16 

64.68 

.0167 

65.19 

.0159 

17 

72.77 

.0171 

72.29 

.0228 

18 

77.06 


77.06 

.0075 

Filter 

— 

— 

32.33 

.316 
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APPENDIX B 


FLUTTER SUPPRESSION SYSTEM IMPLEMENTATION 


This appendix contains data pertaining to the implementation of the redun- 
dant, digital FSS for the B-52E aeroelastic wing tunnel model. 

Detailed circuit diagrams of the analog voter and the analog voter DC power 
supplies are given on Figures B-1 and B-2. 

Listings of the software implementing the full FSS filter and failure 
detection and timing logic are presented on pages B-4 through B-19. A 
listing of the program designed to isolate failed components in the FSS 
is given on pages B-20 through B-22. 

An alphabetic list of assembly language instructions for the HP2100 mini- 
computer is given on pages B-23 through B-25, 

The final analog computer patching diagram is presented on Figure B-3, 
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FIGURE B-1 - ANALOG VOTER CIRCUIT DIAGRAM 
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FIGURE B-2 - ANALOG VOTER POUR SUPPLIES CIRCUIT DIAGRAM 







FLUTTER SUPPRESSION SYSTEM COMPUTER PROGRAM 


0001 ASHB»A,B»T»L 

0002*THIS PR06RAM IS FOR DIGITAL FLUTTER SUPPRESSION SYS FOR B-52 MODEL 
0003AFILE NAME IS " JFIL8 ' CREATED BY JR MATTHEW 8/03/78 
0004*BILINEAR TRANSFORMATION METHOD IS USED TO IMPLEMENT THIS FILTER 
0005*HP 5610A A TO 0 CONVRTR AND HP 6940A MULTIPROG, ARE FOR USEDFOR DATA 10 
0006* SAMPLING PERIOD *.002 MILLISECONDS 


0007*< 






0006* 

0009* 

INITIALIZE 

THE HP 1/0 DEVICES 

0010 

0011 

00100 

00100 

002400 

ORG 

CLA 

1000 

INITIALIZE MEM LOC 12B TO ZERO(NOP) 

0012 

00101 

070012 

3TA 

12B 

0013 

00102 

070011 

3TA 

118 

INITIALIZE MEM LOC IIB 

0014 

00103 

061165 

LO* 

INIT 

LOAD INT REG'A* CONTROL WORD 1701408 

0015 

00104 

102612 

OTA 

123 


0016 

00105 

103712 

STC 

12B,C 


0017*< 



*** «* 

**********4 

TO 1,0 MILLISEC clock PULSES 

0018*INITIALIZE TIME 
0019* 

BASE 

GENERATOR 

0020 

00106 

102100 

STF 

0 


0021 

00107 

061170 

LDA 

Cw 


0022 

00140 

102610 

OTA 

lOB 


0023 

00111 

061173 

LDA 

USB 


0024 

00112 

070010 

STA 

lOR 


0025 

00113 

061075 

LOA 

CNT 


0026 

00114 

070020 

STA 

20B 


0027*^ 


r******** 

***** 

**********^ 

PROGRAM LOOP START AT THIS POINT 

0026 

0029* 

00115 

015066 

START 

JSB RCNTR 

0030*^ 



« *** * 

**********^ 


0031* 

0032* 

0033* 

0034* 

0035 

0036 

S REGISTER FLASH ROUTINE 
FLASHES AT 2HZ 

00116 000000 NOP 

00117 002400 CLA 


0037 

00120 

065303 

LOB 

AFF 

»0»s A CHANNEL FINISHED FLAG 

0038 

00121 

006020 

SSB 


FINISHEOs-1 »NOT FINISHED*© 

0039 

00122 

024125 

JMP 

*♦3 


0040 

00123 

035310 

ISZ 

EPCNT 


0041 

0042 

00124 

00125 

031312 

000000 

lOR 

NOP 

AFFER 

SET BIT 2 IN FLASH I^ORD 

0043 

00126 

065304 

LDB 

BFF 

*8'» B channel FINISHED FLAG 

0044 

00127 

006020 

SSB 


FINISHED*-! , NOT FINISHED*© 

0045 

00130 

024133 

JMP 

*♦3 


0046 

00131 

035310 

ISZ 

ERCNT 


0047 

0046 

00132 

00133 

031313 

000000 

lOR 

NOP 

BFFER 

SET BIT 5 IN FLASH wORO 

0049 

00134 

065306 

LOB 

ACF 

*B*« A CHANNEL COMPARISON FLAG 

0050 

00135 

006020 

SSB 


NO ERR0R*-1 ,ERROR*0 

0051 

00136 

024141 

JMP 

*♦3 


0052 

00137 

035310 

ISZ 

ERCNT 


0053 

0054 

00140 

00141 

031314 

000000 

lOR 

NOP 

ACFER 

SET BIT 6 IN FLASH WORD 

0055 

00142 

065307 

LDB 

BCF 

*B's B CHANNEL COMPARISON FLAG 

0056 

00143 

006020 

SSB 


NO ERROR*-lfCRROR*0 
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FLUTTER 

SUPPRESSION 

SYSTEM COMPUTER PROGRAM (CONTINUED) 

005T 

00144 

024147 

JMF 

*♦3 



0056 

00145 

035310 

ISZ 

ERCNT 



0059 

00146 

031315 

lOR 

BCFER 


SET BIT 11 IN FLASH WORD 

0060 

00147 

000000 

NOP 




0061 

00150 

065305 

LOB 

CFF 


*B'B COMPARISON FINISHED FLAG 

0062 

00151 

006020 

SSB 



FINISHE0«*1»N0T FINISHED«0 

0063 

00152 

024155 

JMP 

*♦3 



0064 

00153 

035310 

ISZ 

ERCNT 



0065 

00154 

031316 

lOR 

CFFER 


SET BIT 14 IN FLASH WORD 

0066 

00155 

000000 

NOP 




0067* 







0066* 

CHECK 

FOR ERRORS AND DO FLASH 

OUTPUT 

0069* 







0070 

00156 

035311 

ISZ 

ERCLK 


INCREMENT ERROR CLOCK SKIP IF « 0 

0071 

00157 

024176 

JMP 

CLEAR 


NOT TO CHECK FOR ERRORS 

0072 

00160 

000000 

NOP 



TIME TO CHECK FOR ERRORS 

0073 

00161 

065310 

LOB 

ERCNT 


»e'« number of ERRORS/SAMPLE PERIOO 

0074 

00162 

006003 

SZB, 

RSS 


ANY ERRORS????? 

0075 

00163 

024170 

JMP 

GOOD 


NO 

0076 

00164 

106501 

Lie 

IB 


YES-GET FLASH wORO FROM s REGISTER 

0077 

00165 

006021 

SSB# 

RSS 


WAS FLASH WORD COMPLEMENTED LAST TIME 

0076 

00166 

003000 

CMA 



no-complement flash word 

0079 

00167 

024171 

JMP 

FLASH 


YES-GO TO OUTPUT ROUTINE 

0060 

00170 

002400 

GOOD 

CLA 


CLEAR FLASH WORD 

0061 

00171 

102601 

FLASH 

1 OTA ! 

LB 

PUT FLASH WORD INTO S REGISTER 

0082 

00172 

065317 

LOB 

FLCNT 


-PCYCLES/FLASH 

0083 

00173 

075311 

STB 

ERCLK 


RESET ERROR CLOCK 

0084 

00174 

006400 

CLB 




0085 

00175 

075310 

STB 

ERCNT 


SET ERROR COUNT TO ZERO 

0086 

00176 

006400 

CLEAR 

t CLB 



0087 

00177 

075303 

STB 

AFF 


CLEAR A channel FINISHED FLAG 

0088 

00200 

075304 

STB 

BFF 


CLEAR 0 channel FINISHEO FLAG 

0089 

00201 

075306 

STB 

ACF 


CLEAR A CHANNEL COMPARISON FLAG 

0090 

00202 

075307 

STB 

BCF 


CLEAR B channel COMPARISON FLAG 

0091 

00203 

075305 

STB 

CFF 


clear channel COMPARISON FLAG 

0092 

00204 

000000 

NOP 




0095*< 







0094*CHANNEL A F3S 

LOOP 




0095*t 



«****ri 

k***A*i 

*****1 


0096*READ INPUT FROM CH ( 

), INPU 

PORT 

llB 

0097* 







0098 

00205 

061157 

LOA 

CHA 

SELECT CHANNEL 0 OF A/0 DEVICE (I/O PORT UR) 

0099 

00206 

102611 

OTA 

IIB 

OUTPUT 'A* REG TO I/O PORT IIB 

0100 

00207 

103711 

STC 

11B,C 


DEVICE command TO I/O llB 

0101 

00210 

102311 

SFS 

1 IB 



0102 

00211 

024210 

JMP 

*-l 



0103 

00212 

102511 

LIA 

IIB 

READ I/O BUFFER CONTENT IN *A' REG 

0104 

00213 

011166 

AND 

MASK 


ZERO OUT/ CH.IO ON BITS 0 THRU 5 

0105 

00214 

071145 

STA 

INPTA 



0106*1 




****** 



0107*CHANNEL A FSS 

FILTER IMPLEMENTATION 

0108*1 







0109* 

STAGE All LAG (10/(6*10)] 


0110* 







0111 

00215 

000000 

NOP 




0112 

00216 

061126 

LOA 

YlAll 


YlAll IN 'A'^YIN-DTm 
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FLUTTER SUPPRESSION SYSTEM COMPUTER PROGRAM (CONTINUED) 


0113 

00217 

100200 

MPY 

KlAll 

♦B'«Y(N-nT*KlAll (2**-l) 



00220 

001204 





0114 

00221 

100021 

ASL 

1 

*e*»y (N-1)*K1A11 


Oils 

00222 

045127 

AOB 

KXAll 

*0'a Y(N-1)*K1A11*X(M.1)T*K2A11 

0116 

00223 

075144 

STB 

TEMP 



0117 

00224 

061145 

LOA 

INPTA 

*A*«X(NT),,, 'INPTA* 

11, 

0110 

00225 

100200 

MPY 

K2A11 

*8'*X(NT)*K2A11 (2**-l) 



00226 

001205 





0119 

00227 

100021 

ASL 

1 

*B'«X(NT)*K2A11 


0120 

00230 

075127 

STB 

KXAll 

KXA11»X(N-1)T*K2A11 FOR NEXT ITERATION 

0121 

00231 

045144 

ADB 

TEMP 

*B**(X (NT)*X (N-1) T*K2A1 1 * 

Y(N-1)T*K1 All 

0122 

00232 

075126 

STB 

YlAll 

Y1A11*Y(NT),, OUTPUT OF THIS 

STAGE 

0123 

00233 

000000 

NOP 




ui2a*i 

»*****# 






0125* 

STAGE A12 LAG (10/(8*10)) 



0126* 







0127 

00230 

000000 

NOP 




0128 

00235 

061130 

LOA 

Y1A12 

Y1A12 IN 'A»*Y (N-l)Tf » 


0129 

00236 

100200 

MPY 

K1A12 

'e**Y(N-l)T*KlA12(2**-n 



00237 

001206 





0130 

00240 

100021 

ASL 

1 

•H'sY(N-1)*k1A12 


0131 

002<M 

045151 

AOB 

KXA12 

'B's Y(N-l)*KlA12*X(N-l)T*i<2A12 

0132 

00202 

075144 

STB 

TEMP 



0133 

00203 

061126 

LDA 

YlAll 

*A**X(NT)=Oi)TPUT from stage 

All 

0134 

00200 

100200 

MPY 

K2A12 

'B*sX(nt)*K2A12(2**-D 



00205 

001207 





0135 

00246 

100021 

ASL 

1 

*H'aX(NT)*<2A12 


0136 

00247 

075131 

STB 

KXA12 

KXA12sX(N-l)T*K2A12 FOR NEXT ITERATION 

0137 

00250 

045144 

AOB 

TEMP 

*()*«(X(NT)*X(N-1) T*t<2A12 + 

Y(N-1)T*K1M2 

0138 

00251 

075130 

STB 

Y1A12 

YlA12rY(.^T) ,, OUTPUT OF THIS 

STAGE 

0139 

00252 

000000 

NOP 





0140****************************** 

0141* STAGE A21 lag 150/(8*150) 
0142* 


0143 

00253 

061132 

LOA 

Y1A21 

* A 'sY (N-l ) 

0144 

00254 

100200 

MPY 

KU21 

'B»aY(N-l)*KlA21(2**-l) 

0145 

00255 

00256 

001210 

100021 

ASL 

1 

'B'aY(N-l)*KlA21 

0146 

00257 

045133 

ADB 

KXA21 

»8*aY(W-n*KlA21*X(N-l)*n2A21 

0147 

00260 

075144 

STB 

TEMP 


0148 

00261 

061145 

LOA 

INPTA 

'A'a X(NT)m|FSS channel A INPUT 

0149 

00262 

100200 

MPY 

K2A21 

'H'xX(NT)*K2A21 (2**»1) 

0150 

00263 

00264 

001211 

100021 

ASL 

1 

*B* X(nT)*«2A21 

0151 

00265 

075133 

STB 

KXA21 

KXA21XX (N-n *K2A21 FOR NEXT ITERATION 

0152 

00266 

045144 

AOB 

TEMP 

♦B'a(X(NT)*X(N«l)T)*K2A21 *V(N- 1 )T*K 1 A 21 

0153 

00267 

075132 

STB 

Y1A21 

YlA21aY(NT) ,, OUTPUT OF THIS STAGE 

0154 

0155*1 

00270 

******* 

000000 

NOP 

r*** ^ 



0156* 

STAGE 

: A22 LAG 

; iso/(s*i5o) 


0157* 

0158 

00271 

061134 

LDA 

Y1A22 

*A*aY(N*l) 

0159 

00272 

100200 

MPY 

K1A22 

'8'aY(N-l)*KlA22(2**-l) 

0160 

00273 

00274 

001212 

100021 

ASL 

1 

*B'aY(N-n*KlA22 

0161 

00275 

045135 

AOB 

KXA22 

*e'aY(N-n*KlA22*X(N*l)*K2A22 
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FLUTTER SUPPRESSION SYSTEM COMPUTER PROGRAM (CONTINUED) 


0162 

00276 

075144 

ST8 TEMP 

•A'« X(NT)«OUTPUT FROM STAGE A21 

0163 

00277 

061132 

LDA Y1A21 

0164 

00300 

00301 

100200 

001213 

MPY K2A22 

*0**X(NT)*K2A22(2**-l) 

0165 

00302 

100021 

ASL 1 

*8* X(NT)*K2A22 

0166 

00303 

075135 

STB KXA22 

KXA22«X(N-l)*K2A22 FOR NEXT ITERATION 

0167 

00304 

045144 

AOB TEMP 

•0*B(X(NT)*X(N»nT)*K2A22 ♦Y(N-1)T*K1A22 

0166 

00305 

075134 

STB Y1A22 

Y1A22*Y(NT),, OUTPUT OF THIS STAGE 

0169 

00306 

000000 

NOP 


0170** 

lr**«**« 



0171* 

0172* 

STAGE 

A23 LAG 150/(3*150) 

*A*«Y(N«n 

0173 

00307 

061136 

LDA YIA23 

0174 

00310 

00311 

100200 

001214 

MPY K1A23 

•8»»Y(N-1)*KIA23(2**-1) 

0175 

00312 

100021 

ASL 1 

•8'«Y(N-1)*K1A23 

0176 

00313 

045137 

AOB KXA23 

»B'»V(N-l)*KlA23*X(N.n*K2A23 

0177 

00314 

075144 

3T0 TEMP 

«A's X(NT)»0UTPUT from STAGE A 3 ? 

0176 

00315 

061134 

LDA Y1A22 

0179 

00316 

00317 

100200 

001215 

MPY X2A23 

*H'xx (NT) *K2A23(2**-1 ) 

0160 

00320 

100021 

ASL 1 

'0' X(NT)*K2A23 

0161 

00321 

075137 

STB KXA23 

KXA23ax(N-n»t<2A23 FOR NEXT ITERATION 

0162 

00322 

045144 

AD8 TEMP 

*e*a(X(NT)*X(N-l)T)*K2A23 *Y(M-1)T*K1A23 

0163 

00323 

075136 

STB Y1A23 

Y1A23*Y (NT) f fOUTPUT OF THIS STAGE 

0164 

00324 

000000 

NOP 


0185* 





0186** STAGE A31 F IRST-QROER/SECOND 

-OHOEH 

0187** 32,1325(S*5b,8in / (S2*50S*15625) 

0188** 

0169 

00325 

000000 

NOP 

*A'*Y(N-2) 

0190 

00326 

061141 

LOA Y2A31 

0141 

00327 

00330 

100200 

001217 

MPY K2A31 


0192 

00331 

100021 

ASL 1 

»B*=Y(N-2)*K2A31 

0193 

00332 

075150 

STB TEMPI 

»A*=Y(N-1) 

0194 

00333 

061140 

LOA Y1A31 

0195 

00334 

071141 

STA Y2A31 

PUT Y(N-l) INTO Y(N-2) FOP NEXT IT 

0196 

00335 

00336 

100200 

001216 

MPY K1A31 

»6*sY(N-n*FlA31 (2**-2) 

0197 

00337 

100021 

ASL 1 

»P'*Y(N-1)*K1A31(2**-1) 

0196 

00340 

075151 

3TB TEMP2 

»A*aX(N-2) 

0199 

00341 

061143 

LOA X3A31 

0200 

00342 

00343 

100200 

001222 

MPY K5A31 


0201 

00344 

100021 

ASL 1 

'B'xX (N4r 2) wKSASl 

0202 

00345 

075152 

STB TEMP3 

*A'sX(N-1) 

0203 

00346 

061142 

LOA X2A31 

0204 

00347 

071143 

STA X3A31 

PUT X(N-l) INTO X(N-2) FOP NEXT IT 

0205 

00350 

00351 

100200 

001221 

MPY K4A31 


0206 

00352 

100021 

ASL i 

'F««X(N-1)*K4A31 

0207 

00353 

075153 

STB TEMP4 

•A»«X(N),,,,FSS channel a input 

0208 

00354 

061145 

LDA INPTA 

0209 

00355 

001121 

ARSr ARS 

*A*«X(N)/8 PREVENT OVERFLOW AT 2 

0210 

00356 

001100 

ARS 
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FLUTTER SUPPRESSION SYSTEM COMPUTER PROGRAM (CONTINUED) 


0211 

00357 

071142 

STA X2A31 

0212 

00360 

100200 

MPY K3A31 


00361 

001220 


0213 

00362 

100021 

ASL 1 

021tt 

00363 

045152 

AOB TEMP3 

0215 

00364 

000000 

NOP 

0216 

00365 

045153 

AOB TEMP4 

0217 

00366 

000000 

NOP 

0216 

00367 

045150 

ADB TEMPI 

0219 

00370 

000000 

NOP 

0220 

00371 

005100 

8R5 1 

0221 

00372 

045151 

AOB TEMP2 

0222 

00373 

000000 

NOP 

0223 

00374 

005000 

BLS 1 

0224 

00375 

075140 

STB Y1A31 

0225 

00376 

000000 

NOP 

0226«< 



************* 

0227* 

SUMMATION OF 

FILTER OUTPUT 

0228* 




0229 

00377 

061126 

LOA YlAll 

0230 

ooaoo 

100200 

MPY KOAll 


ooaoi 

001224 


0231 

ooao 2 

100021 

A3L 1 

0232 

00403 

075150 

STB TEMPI 

0233 

oo4oa 

061130 

LD* Y1A12 

0234 

00405 

100200 

MPY KOAIP 


00406 

001223 


0235 

00407 

100021 

ASL 1 

0236 

00410 

075151 

STB TEMP2 

0237 

0041 1 

061132 

LDA YIA21 

0236 

00412 

100200 

MPY K0A21 


00413 

001227 


0239 

00414 

100021 

ASL 1 

0240 

00415 

075152 

STB TEMP3 

0241 

00416 

061134 

LO* Y1A22 

0242 

00417 

100200 

MPY KUA22 


00420 

001226 


0243 

00421 

100021 

ASL 1 

0244 

00422 

075153 

STB TEMP4 

0245 

00423 

061136 

LO* Yl*23 

0246 

00424 

100200 

MPY K0A23 


00425 

001225 


0247 

00426 

100021 

ASL 1 

0248 

00427 

075154 

STB TEMPS 

0249 

00430 

061140 

LOA YIA31 

0250 

00431 

100200 

MPY K0A31 


00432 

001230 


0251 

00433 

100021 

ASL 1 

0252 

00434 

045154 

AOB TEMP5 

0253 

00435 

045153 

ADB TEMP4 

0254 

00436 

045152 

ADB TEMP3 

0255 

00437 

045151 

ADR TEMP2 

0256 

00440 

045150 

AOB TEMPI 

0257 

00441 

103201 

SOC C 

0256 

00442 

015106 

JSB OVFLO 

0259 

00443 

075155 

STB OUT* 


PUT X(N) INTO X(N*-1) FOR NEXT ITER 

*R'»X(N)*K3A31 

*»•« " ♦ X(N-2)*K5A31 

»B*« " ♦ X(N-l)*K«*3l 

*B»a " ♦ Y(N-2)*K2A31 

♦b'« "(2**-1) 

*P'» ■ ♦Y(N-1)*K1A31 

•B'« Y (N) ' t know last 

Y1A3I3Y(NT) ,, OUTPUT OF THIS STAGE 


•A'sOUTPUT from stage All 

'B's All*KOAll 
'A'sOUTPL'T FROM STAGE A12 

*«•= A12*K0A12 
»4»a0uTPUT FROf^ STAGE A21 

•P'a A21*KOA2l 
'a'sOUTPUT from stage A22 

*P('= 422*K0A22 
•A'aOUTPUT FROM STAGE *23 

'S'a A 23 *«OA 23 
'A'aOUTPUT FROM STAGE A 31 

'8»a A3l*K0A31 
•8'* " ♦A23*K0*23 
♦B'a ■ +A22AK0A22 
«n*a * ♦A21*K0A21 
•B's " ♦A12*K0A12 
•B'a " ♦All+KOAll 

SATURATE IF OVERFLOW OCCURRS 
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FLUHER SUPPRESSION SYSTEM COMPUTER PROGRAM (CONTINUED) 


0260 

00444 000000 

NOP 


0261** 


**** 


0262*CHANNEL A OUTPUT 

ROUTINE 

0263* 




0264*APPLT PREDICTOR 

algorithm 

0265* 




0266* 

Y (N)«KIP1*X(N)*K2P1*X (N-1 

0267* 




0266 

00445 000000 

NOP 


0269 

00446 061147 

LOA 

INPA2 

0270 

00447 100200 

MPY 

K3PA1 


00450 001203 



0271 

00451 100021 

ASL 

\ 

0272 

00452 075150 

3TR 

TEMPI 

0273 

00453 061146 

LOA 

INPAl 

0274 

00454 071147 

STA 

IMPA2 

0275 

00455 100200 

MPY 

K2PA1 


00456 001202 



0276 

00457 100021 

ASL 

1 

0277 

00460 075144 

STB 

TEMP 

0278 

00461 061155 

LDA 

OUTA 

0279 

00462 071146 

STA 

INPAl 

0280 

00463 100200 

MPY 

KlPAl 


00464 001201 



0281 

00465 100021 

ASL 

1 

0282 

00466 045144 

AOB 

TEMP 

0283 

00467 045150 

ADB 

TEMPI 

0284 

00470 100021 

ASL 

1 

0285 

00471 103201 

30C 

C 

0286 

00472 015106 

JS6 

OVFLO 

0287* 




0288*0UTPUT CHANNEL 

A ON 

0/A Cl 

0289* 




0290 

00473 101044 

LSR 

a 

0291 

00474 102312 

SfS 

120 

0292 

00475 024474 

JMP 

*-l 

0293 

00476 106612 

0T8 

12B 

0294 

00477 103712 

8TC 

12BfC 

0295 

00500 000000 

NOP 


0296* 





TO OUTPUT SIGNAL 


»A*aOUTA(N-2) 


*B**OUTA (N-2)*K3PAl (?**■! ) 


* 

PUT 


!0UTA(N-1) 

OUTA(N-l) 


INTO OUTAlN-2) FOR NEXT I 


'0'*OUTA(N-1)*K2PA1(2**-1) 


* A* 

PUT 


tOUTA(N) 

OUTA(N) 


INTO outa(n-i) for iteraion 


'B**0UT(N)*K1P1 

*«'« " ♦ OUT A (N»l ) *F?PA 1 (2**-l ) 

•R*« " ♦ OUTA (N-?) *K3PAl (2**-n 
*0** PREDICTED OUTPUT .002 MILLISFC 
CHECK FOR OVERFLOW 


AH 


ZERO 


LSR to include CH 008 ID, IN BITS 12 THRU IS 

12 BIT DATA A aeiT CH ID. TO OUTPUT BUFFER I 


0297* 

0298* 

0299 

0300 


SET channel a FINISHED FLAG 


00501 061320 LDA FLAG 

00502 071303 STA AFF 

0301 00503 000000 NOP 

0302*********************************** 
0303WCHANNEL B FSS LOOP 
0304*********************************** 

0305WREAD INPUT FROM CH 1»INPU PORT IIB 

0306* 


0307 

00504 

061160 

LOA 

CHB 

0308 

00505 

102611 

OTA 

IIB 

0309 

00506 

103711 

STC 

11B,C 

0310 

00507 

102311 

SFS 

118 

0311 

00510 

024507 

JMP 

*-l 

0312 

00511 

102511 

LIA 

118 


SELECT CHANNEL I OF A/D DEVICE 
OUTPUT 'A' BEG TO I/O PORT IIB 
DEVICE command TO I/O llB 


(I/O PORT 13B) 


READ I/O BUFFER CONTENT IN *A* REG 
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FLUTTER SUPPRESSION SYSTEM COMPUTER PROGRAM (CONTINUED) 

03t3 00512 Olllfcti AND MASK ZERO OUT/ CH.ID ON BITS 0 THRU 5 

0314 00513 071247 STA INPTB 

0315*************************************************************** 

0316*CHANNEL B F88 FILTER IMPLEMENTATION 


0317*< 






031B* 

STAGE 

Bll LAG 

(10/(S*10)) 


0319* 






0320 

oosia 

000000 

NOP 



0321 

00515 

061231 

LOA 

YlBll 

YlBll IN *A»xV(N-l)T 

0322 

00516 

100200 

MPY 

KlBll 

»H»*Y(N-1)T*K1B11 (2***l) 


00517 

001256 




0323 

00520 

100021 

ASL 

1 

*8*=Y(N«1)*K1B11 

0324 

00521 

045232 

AOB 

KXBl 1 

•B's Y(N-1)*KIB11*X(N-1)T*K2R11 

0325 

00522 

075144 

STB 

TEMP 


0326 

00523 

061247 

LDA 

INPTB 

♦A'*X(NT),,,»INPTB» 

0327 

00524 

100200 

MPY 

K2B11 

•B*»X(NT)*K2Bll(2**-l) 


00525 

001257 




0320 

00526 

100021 

ASL 

1 

*8*sX(NTl*K2Bll 

0329 

00527 

075232 

STB 

KXBll 

KXBll*X(N-l)T*K20ll FOR NEXT ITERATION 

0330 

00530 

045144 

AOB 

TEMP 

'B'«(X(NT)*XfN-nT*K2Bll + Y(N-1)T*K1B1 

0331 

00531 

075231 

STB 

YlBll 

Y1B11*Y(NT),, OUTPUT OF THIS STAGE 

0332 

00532 

000000 

NOP 



0333*^ 

»****«« 


r * ★ ★ i 

^^*********** 


0334* 

STAGE 

: 812 LAG (10/(S*10)) 


0335* 






0336 

00533 

000000 

NOP 



0337 

00534 

061233 

LDA 

Y1B12 

Y1812 IN *A'aY(N-l)T 

0330 

00535 

100200 

MPY 

K1P12 

'8'sY(N-nT*KlP12(2**-l) 


00536 

001260 




0339 

00537 

100021 

ASL 

1 

»p'aY(N»l)*KlB12 

0340 

00540 

045234 

AOB 

KXB12 

'B's Y(N-n*KlB12*X(N-.l)T*K2B12 

0341 

00541 

075144 

STB 

TEMP 


0342 

00542 

061231 

LDA 

YlBll 

♦A*aX(NT)aOUTPUT FROM STAGE 811 

0343 

00543 

100200 

MPY 

K2ei2 

•B'aX (NT) *K2fll2(2**-l) 


00544 

001261 




0344 

00545 

100021 

A8L 

1 

«8'xX(NT)*K2B12 

0345 

00546 

075234 

STB 

KXB12 

KXB12*X(N-nT*K2B12 FOR NEXT ITERATION 

0346 

00547 

045144 

AOB 

TEMP 

*B'*(X(NT)+X(N*1)T*K2P12 ♦ Y ( N- 1 ) T*K 1 B 1 2 

0347 

00550 

075233 

STB 

Y1B12 

YlB12eY (NT) ,, OUTPUT OF THIS STAGE 

0348 

00551 

000000 

NOP 



0349*< 


r«*****«*^ 


^***********<* 


0350* 

STAGE B21 LAG l50/fS+l50) 


0351* 






0352 

00552 

061235 

LDA 

YIB21 

*A'aV(N-l) 

0353 

00553 

100200 

MPY 

K1B21 

*fl*aY(N-l)*KlB2l (2**-l) 


00554 

001262 




0354 

00555 

100021 

ASL 

1 

*B»aV(N-l)*KlB2l 

0355 

00556 

045236 

AOB 

KXB21 

*8*aY(N-l)*K1821*X(N-l)*K2B21 

0356 

00557 

075144 

STB 

TEMP 


0357 

00560 

061247 

LOA 

INPTB 

•A»x X(NT)M|F8S channel B INPUT 

0358 

00561 

100200 

MPY 

K2B21 

*B»«X(NT)*K2B21 (2***1) 


00562 

001263 




0359 

00563 

100021 

ASL 

1 

'B' X(NT)*K2B21 

0360 

00564 

075236 

STB 

KXB21 

KXB21aX(N«l)*K2B21 FOR NEXT ITERATION 

0361 

00565 

045144 

AOB 

TEMP 

*B»«(X(NT)*X(N-1)T)*K2B21 ♦ V (N-1 ) T*K 1B2 1 

0362 

00566 

075235 

STB 

Y1B21 

Y1B21»Y(NT)»»0UTPUT OF THIS STAGE 


B-IO 



FLUTTER SUPPRESSION SYSTEM COMPUTER PROGRAM (CONTINUED) 


0363 

00567 

000000 

NOP 




****** 


i***i 



0365* 

STAGE 

: 822 LAG 

1 150/(8*150) 


0366* 






0367 

00570 

061237 

LDA 

riB22 

•A»«Y (N*1 ) 

0369 

00571 

100200 

MPY 

KIB22 

♦B»«Y(N-1)*K1B22(2**-1) 


00572 

001264 




0369 

00573 

100021 

ASL 

1 

•B»*Y(N-1)*KIB22 

0370 

00574 

045240 

A08 

KX622 

•e»«Y (N-1)*K1B22*X(N«1 )*k2822 

0371 

00575 

075144 

ST8 

TEMP 


0372 

00576 

061235 

LOA 

Y1B21 

X(NT)*0UTPUT from stage P21 

0373 

00577 

100200 

MPY 

K2B22 

*8'«X(MT)*«2B22(2**-1) 


00600 

001265 




0374 

00601 

100021 

ASL 

1 

'8' X(NT)*K2B22 

0375 

00602 

075240 

STB 

KX622 

KXB22*X(N-1)*K2B22 FOR NEXT ITERATION 

0376 

00b03 

045144 

ADB 

TEMP 

*e*«(X(NT)*X(N-l)T)*K2B22 *Y (N-l ) T*K 1B22 

0377 

00604 

075237 

ST0 

Y1B22 

Y1«22«Y (NT) ,, OUTPUT OF THIS STAGE 

0378 

00605 

000000 

NOP 



0379** 



r * * * 1 

1*********** 


0380* 

STAGE 

: B23 LAG 

; 150/(8*150) 


0381* 






0382 

00606 

061241 

LOA 

YIB23 

•A'sY(N-I) 

0383 

00607 

100200 

MPY 

K1023 

'0’aY(N-l)*KlB23(2**-l) 


00610 

001266 




0384 

0061 1 

100021 

ASL 

1 

•B'aY(N-n*KlB23 

0385 

00612 

045242 

A08 

KXB23 

'0*sY(N-n*KlB23*X(N-n*K2B23 

0386 

00613 

075144 

STS 

TEMP 


0387 

00614 

061237 

LOA 

YIB22 

X(NT)aOUTPUT FROM STAGE B22 

0388 

00615 

100200 

MPY 

K2R23 

'H'aX(NT)*K2B23(2**-D 


00616 

001267 




0389 

00617 

100021 

ASL 

1 

*B' X(NT)*K2B23 

0390 

00620 

075242 

ST8 

KXB23 

KX823aX(N-l)*K2023 FOR NEXT ITERATION 

0391 

00621 

045144 

AOB 

TEMP 

*H'a(X(NT)*X(N-l)T)*K2B23 *Y (N-l ) T*K1B23 

0392 

00622 

075241 

STB 

Y1B23 

Y 1P23SY (NT) , fOUTPUT OF THIS STAGE 

0393 

00623 

000000 

MOP 



0394** 



r * * * ] 



0395** 

r STAGE B3l FIRST^I 

:iRDER/SECOND 

-ORDER 

0396** 32,1325(S*56.I 

B11)/(S2*50S*15625) 

0397** 





0398 

00624 

000000 

NOP 



0399 

00625 

061244 

LDA 

Y2831 

•a»3Y(N-2) 

0400 

00626 

100200 

MPY 

K2B31 



00627 

001271 




0401 

00630 

100021 

ASL 

1 

•e*ay(N-2)*K2B31 

0402 

00631 

075150 

ST8 

TEMPI 


0403 

00632 

061243 

LOA 

Y1B31 

«A'aY(N-l) 

0404 

00633 

071244 

STA 

Y2B31 

PUT Y(N-l) INTO Y(N-2) FOR NEXT IT 

oaos 

00634 

100200 

MPY 

K1B31 

♦B'aY (N-l) *K1B31 (2**-2) 


00635 

001270 




0406 

00636 

100021 

ASL 

1 

»B»aY(N-n*KlB31 (2**-l) 

0407 

00637 

075151 

ST8 

TEMP2 


0408 

00640 

061246 

LOA 

X3B31 

♦A'ax (N-2) 

0409 

00641 

100200 

MPY 

K5B31 



00642 

001274 




0410 

00643 

100021 

ASL 

1 

»B'aX(N-2)*K5B31 

0411 

00644 

075152 

STB 

TEMP3 
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FLUTTER SUPPRESSION SYSTEM COMPUTER PROGRAM (CONTINUED) 


0412 

00645 

061245 

LOA 

X2B31 

0413 

00646 

071246 

ST4 

X3B31 

04ia 

00647 

00650 

100200 

001273 

MPY 

K4B31 

0415 

00651 

100021 

ASL 

1 

04U> 

00652 

075153 

STB 

TEMP4 

0417 

00653 

061247 

LOA 

INPTB 

0418 

00654 

001121 

ARS, 

ARS 

0419 

00655 

001100 

AR9 


0420 

00656 

071245 

STA 

X2831 

0421 

00657 

00660 

100200 

001272 

MPY 

K3B31 

0422 

00661 

100021 

ASL 

1 

0423 

00662 

045152 

AOB 

TEMP3 

0424 

00663 

000000 

NOP 


0425 

00664 

045153 

AOB 

TEFP4 

0426 

00665 

000000 

NOP 


0427 

00666 

045150 

AOB 

TEMPI 

0426 

00667 

000000 

NOP 


0429 

00670 

005100 

BRS 

1 

0430 

00671 

045151 

AOB 

TEMP2 

0431 

00672 

000000 

NOP 


0432 

00673 

005000 

8LS 

1 

0433 

00674 

075243 

3T8 

Y1B31 

0434 

00675 

000000 

NOP 


0435*1 




^* **** * 

0436* 

0437* 

SUMMATION OF 

filter outpu 

0438 

00676 

061231 

LOA 

YlBll 

0439 

00677 

00700 

100200 

001276 

MPY 

KOBl 1 

0440 

00701 

100021 

ASL 

1 

0441 

00702 

075150 

STB 

TEMPI 

0442 

00703 

061233 

LDA 

Y1B12 

0443 

00704 

00705 

100200 

001275 

MPY 

K0812 

0444 

00706 

100021 

ASL 

1 

0445 

00707 

075151 

STB 

TEMP2 

0446 

00710 

061235 

LDA 

V1821 

0447 

00711 

00712 

100200 

001301 

MPY 

K0B21 

0448 

00713 

100021 

ASL 

1 

0449 

00714 

075152 

STB 

TEMP3 

0450 

00715 

061237 

LDA 

V1B22 

0451 

00716 

00717 

100200 

001300 

MPY 

K0B22 

0452 

00720 

100021 

ASL 

1 

0453 

00721 

075153 

STB 

TEMP4 

0454 

00722 

061241 

LOA 

VIB2S 

0455 

00723 

00724 

100200 

001277 

MPY 

K0B23 

0456 

00725 

100021 

ASL 

1 

0457 

00726 

075154 

STB 

TEMP5 

0458 

00727 

061243 

LOA 

Y1031 

0459 

00730 

00731 

100200 

001302 

MPY 

K0B31 


• A'«x(N*n 

PUT X(N»l) INTO X(N-2) FOR NEXT IT 


•B'aX(N*n*Kae31 

• A'sX(N) r f f »FSS CHANNEL B INPUT 

*A*sX(N)/8 prevent overflow *T 2 
PUT X(N) INTO X(N»1) FOR NEXT ITER 


•B'*X(N)*K3051 

*B*« ■ ♦ X(N-2)*K5B31 

" ♦ X(N-l)*KaP31 

•0»« " ♦ V(N*2)*K263l 

'B'« "(2***n 

*B*» " ♦V(N-1)*K1B31 C2**-n 

*B*« t(N)»»*»—»*»-DON'T KNOW LAST 
VIP31=V(NT) ,, OUTPUT Of- THIS STAGE 


»a**output from stage en 

•H'a Bll*KOBtl 
•A'aOUTPUT FROM STAGE 812 

'B's 012*KO812 
'A'rOUTPUT FROM STAGE B21 

'«'* B2l*K0P2l 
•A'aOUTPuT FROM STAGE P22 

*B'a B22AK0B22 
•A'aOUTPUT FROM STAGE B23 

'B*« B2S*K0623 
*A»«0UTPUT from stage B3l 
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FLUTTER SUPPRESSION SYSTEM COMPUTER PROGRAM (CONTINUED) 


0460 

00732 

100021 

ASL 

1 

0461 

00733 

045154 

AOB 

TEMP5 

0462 

00734 

045153 

ADB 

TEMP4 

0463 

00735 

045152 

AOB 

TEMP3 

0464 

00736 

045151 

AOB 

TEMP2 

0465 

00737 

045150 

ADB 

TEMPI 

0466 

00740 

103201 

soc 

C 

0467 

00741 

015106 

JS8 

OVFLO 

0468 

00742 

075252 

STB 

0UT8 

0469 

00743 

000000 

NOP 


0470* 




k***** 

0471*CHANNEL B OUTPUT ROUTINE 


»8'« 

•B*a 


eSl*KOBSl 
" tBaj-mOBZS 

• ♦B22*K0B22 
" ♦B21*K0B21 

• ♦B12*K08!2 

• ♦BU+KOBll 

SATURATE IF OVERFLOW OCCURRS 


0472* 

0a73*APPLT 

0474* 

0475* 

0476* 


PREDICTOR ALGORITHM TO OUTPUT SIGNAL 
Y(N)aKlPl*X(N)*KaPl*X(N-l)*K3Pl*X(N"2) 


0477 

00744 

000000 

NOP 



047B 

00745 

061251 

LOA 

INPB2 

'A*a0UTB(N*2) 

0479 

00746 

100200 

MPY 

K3PB1 



00747 

001255 




04B0 

00750 

100021 

ASL 

1 

*0'sOUTB(N-2)*K3P81 (2**-n 

04B1 

00751 

075150 

STB 

TEMPI 


04B2 

00752 

061250 

LOA 

INPBl 

*A'*OuTB(N"l) 

0483 

00753 

071251 

STA 

INPB2 

PUT OUTB(N-n INTO 0UT8(N-< 

04B4 

00754 

100200 

MPY 

K2PB1 



00755 

001254 




04B5 

00756 

100021 

ASL 

1 

*8’*0UTB(N-l)*K2PBl C2**-!) 

0466 

00757 

075144 

STB 

TEMP 


0487 

00760 

061252 

LDA 

OUTB 

* A'aOUT8(N) 

0486 

00761 

071250 

STA 

INPBl 

PUT OUTfl(Nl INTO OUTB(N-l) 

0489 

00762 

100200 

MPY 

KlPBl 



00763 

001253 




0490 

00764 

100021 

ASL 

1 

*H*aOUTb(N)*KlPRl (2**-l) 

0491 

00765 

045144 

AOB 

TEMP 

" ♦ OUT0(N»l)*K2P8l (2 

0492 

00766 

045150 

AOB 

TEMPI 

»8»s " ♦ 0UTB(N«2)*K3PBI (2 

0493 

00767 

100021 

ASL 

1 

'B'a PREDICTED OUTPUT ,00( 

0494 

00770 

103201 

SOC 

C 

CHECK FOR OVERFLOW 

0495 

00771 

015106 

JSB 

OVFLO 


0496*^ 

k*** * *4 


^* * * *1 



0497*0UTPUT 

CHANNEL 

6 ON 

D/A channel one 

0498* 






0499 

00772 

002404 

CLA 

» INA 


0500 

00773 

101104 

RPR 

4 ROTATE 

TO INCLUDE CH IB ID, IN BITS 

0501 

00774 

102312 

3FS 

12B 


0502 

00775 

024774 

JMP 

*-i 


0503 

00776 

106612 

0TB 

12B 12 

BIT DATA R 48IT CH ID, TO OU 

0504 

00777 

103712 

STC 

12B,C 


0505 

01000 

000000 

NOP 



0506*' 






0507* 

SET CHANNEL B 

FINISHED FLAG 


0508* 






0509 

01001 

061320 

LOA 

FLAG 


0510 

01002 

071304 

STA 

BFF 


0511 

01003 

000000 

NOP 



0512* 

****** 






-2) FOR NEXT T 


FOR ITERAION 


► - 1 ) 

MILLISEC 


AH 


12 Thru i5 


BUFFER I 
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FLUTTER SUPPRESSION SYSTEM COMPUTER PROGRAM (CONTINUED) 


0513* 0I6ITAL CROSS CHECK OF A S B CHANNELS ON ANOTHER COMPUTER'S 

051A* OUTPUT VERSUS OUTPUT OF ANALOG VOTER 

0515* OATA IS SHIFTED TO THE RIGHT TO ELIMINATE ALL BUT 

0516* THE SIGN BIT AND 2 SIGNIFICANT BITS 

0517* 

0518 OlOOa 061161 LOA CH2 

0519 01005 102611 OTA IIB SELECT CHANNEL 2 OF A/0 

0520 01006 103711 STC 118, C 

0521 01007 102311 3FS 118 

0522 01010 025007 JMP *.i 

0523 01011 102511 LIA 118 'A'a CHANNEL A VOTER OUTPUT 

0529 01012 065162 LOB CH3 

0525 01013 106611 OTB 118 SELECT CHANNEL 3 OF A/D 

0526 01014 103711 STC 118, C 

0527 01015 102311 3FS 118 

0528 01016 025015 JMP *-i 

0529 01017 106511 LIB 118 *8*» CHANNEL A COMPUTER OUTPUT 

0530 01020 015116 JS8 PROCS PROCESS OATA 

0531 01021 075321 3T8 CHAER 

0532 01022 061163 LOA CH4 

0533 01023 102611 OTA 118 SELECT CHANNEL « OF A/0 

0534 01024 103711 STC IIB,C 

0535 01025 102311 SFS 110 

0536 01026 025025 JMP *•! 

0537 01027 102511 LIA 118 *A*» CHANNEL B VOTER OUTPUT 

0538 01030 065164 LOB CHS 

0539 01031 106611 OTB IIB SELECT CHANNEL 5 OF A/D 

0540 01032 103711 STC 118, C 

0541 01033 102311 SFS 118 

0542 01034 025033 JMP 

0543 01035 106511 LIB 118 »9*« CHANNEL A COMPUTER OUTPUT 

0544 01036 015116 JS9 PROCS PROCESS DATA 

0545 01037 075322 STB CHBER 

0546* 

0547* CHECK FOR EXCESSIVE ERROR AND SET APPROPRIATE FLAGS 
0548* 


0549 

01040 

061321 

LD* 

CHAER 

0550 

01041 

002002 

SZ4 


0551 

01042 

025045 

JMP 

*♦3 

0552 

01043 

061320 

LOA 

FLAG 

0553 

01044 

071306 

STA 

ACF 

0554 

01045 

000000 

NOP 


0555 

01046 

061322 

LOA 

CHBER 

0556 

01047 

002002 

SZA 


0557 

01050 

025053 

JMP 

*♦3 

0550 

01051 

061320 

LOA 

FLAG 

0559 

01052 

071307 

STA 

BCF 

0560 

01053 

000000 

NOP 


0561 

01054 

061320 

LOA 

FLAG 


'A*» CHAN A VOTER - CHAN A OUTPUT 
IS ERROR ZERO 
NO 
TES 

SET CHANNEL A »0K* FLAG 
DON'T SET CHANNEL '0K» FLAG 
'A** CHAN a VOTER - CHAN 8 OUTPUT 
IS ERROR ZERO 
NO ' 

TES 

SET channel B *0k' flag 
DON'T SET channel B 'OK' FLAG 


0562 01055 071305 3TA CFF 

0563 01056 000000 NOP 


SET COMPARISON FINISHED FLAG 


0564** DUMMY COUNTER ROUTINE 4,9MICR0S£C0N08 COUNT 

0565 01057 002400 CLA 

0566 01060 071176 STA NUMBR 

0567 01061 035176 ISZ NUMBR SKIP IF NUMBRaO, NUMBRa«32768 

0568 01062 025061 JMP *•! JUMP TO INCREMENT 
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FLUTTER SUPPRESSION SYSTEM COMPUTER PROGRAM (CONTINUED) 


0569 01063 000000 NOP 

0570 01064 10a044 HLT 44B 

0571 01065 025065 END JMP * 

0572* 

0573* RESET THE COUNTER 

0574 01066 000000 RCNTR NOP 

0575 01067 061174 LOA MNUM 

0576 01070 071175 ST* CNTON 

0577 01071 061176 LOA NUMBR 

0578 01072 071156 STA LCNT SAVE OUTPUT OF DUMMY COUNTER 

0579 01073 103710 STC 108, C 

0580 01074 125066 JMP RCNTR, I 

0581* 

0582* INCREMENT THE COUNTER 

0583 01075 001076 CNT OEF CNTR 

0584 01076 000000 CNTR NOP 

0585 01077 106710 CLC lOB 

0586 01100 035175 ISZ CNTDN 

0587 OllOl 025104 JMP **3 

0588 01102 102077 HLT 778 HALT AFTER ONE PASS * CHECK DATA IN 'NUM8R' 

0589 01103 024115 JMP START 

0590 01104. 103710 STC 10B,C 

0591 01105 125076 JMP CNTR, I 

0592* 

0593********************************************************************* 

0594* THE OVERFLOW SUBROUTINE 
0595* THE input IS IN THE B REGISTER 

0596 01106 000000 OVFLO NOP 

0597 01107 006020 SSB 

0598 OHIO 025113 JMP * + 3 

0599 Ollll 065171 LOB PLUS 

0600 01112 025114 JMP *+2 

0601 01113 065172 LOB minus 

0602 01114 103101 CLO 

0603 01115 125106 JMP OVFLO, I 



0606* SUBROUTINE TO SUBTRACT A-REG FROM B-RE6 AND TAKE ABSOLUTE 
0607* VALUE AND SHIFT DATA RIGHT 12 PLACES 


0608* 

0609 01116 000000 PR0C3 NOP 

0610 01117 003004 CMA,INA CHANGE SIGN ON A>RE6 

0611 01120 044000 ADB OB ADO A«REG TO 6-REG 

0612 01121 006020 SSB H-REG NEGATIVE ? 

0613 01122 007004 CMB,INB YES* CHANGE SIGN 

0614 01123 000000 NOP NO 

0615 01124 101055 LSR 13 SHIFT RIGHT 13 TO SEE IF ERROR IS BIG( 

0616 01125 125116 JMP PROCS,I RETURN 

0617* 

0618* FILTER INITIAL CONDITIONS CHANNEL A 
0619* 

0620 01126 000000 YlAll OCT 0 

0621 01127 000000 KXAH OCT 0 

0622* 

0623 01130 000000 Y1A12 OCT 0 

0624 01131 000000 KXA12 OCT 0 
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FLUTTER SUPPRESSION SYSTEM COMPUTER PROGRAM (CONTINUED) 


0635 * 


0626 

01132 

000000 

V1A21 OCT 0 


0627 

01133 

000000 

KXA21 OCT 0 


0628* 

0624 

01134 

000000 

YIA22 OCT 0 


0630 

01135 

000000 

KXA22 OCT 0 


0631* 

0632 

01136 

000000 

Y1A23 OCT 0 


0633 

01137 

000000 

KXA23 OCT 0 


0634* 

0635 

01140 

000000 

Y1A31 OCT 0 


0636 

01141 

000000 

Y2A31 OCT 0 


0637 

01142 

000000 

X2A31 OCT 0 


0638 

01143 

000000 

X3A31 OCT 0 


0639* 

0640 

01144 

000000 

TEMP OCT 0 


0641 

01145 

000000 

INPTA OCT 0 


0642 

01146 

000000 

INPAl OCT 0 


0643 

01147 

000000 

INPA2 OCT 0 


0644 

01150 

000000 

TEMPI OCT 0 


0645 

01151 

000000 

TEMP2 OCT 0 


0646 

01152 

000000 

TEMP3 OCT 0 


0647 

01153 

000000 

TEMP4 OCT 0 


0646 

01154 

OOOOOO 

TEMP5 OCT 0 


0649 

01155 

000000 

OUTA OCT 0 


0650 

01156 

OOOOOO 

LCNT nCT 0 


0651*i 

k****«^ 

^******* 

******************** 


0652* 

CONSTANTS FOt^ I/O AND ,004 CYCLE 

TIME 

0653* 

0654 

01157 

100000 

CHA OCT 100000 


0655 

01160 

100001 

CH6 OCT 100001 


0656 

01161 

10O002 

CH2 OCT 100002 


0657 

01162 

100003 

CH3 OCT 100003 


0658 

01163 

100004 

CH4 OCT 1000O4 


0659 

01164 

100005 

CH5 OCT 100005 


0660 

01165 

170140 

IN IT OCT 170 140 


0661 

01166 

177700 

MASK OCT 177700 


0662 

01167 

olouoo 

OACH OCT 010000 


0663 

01170 

000001 

OCT 1 


0664 

01171 

077777 

PLUS OCT 77777 


0665 

01172 

100000 

MINUS OCT 100000 


0666 

01173 

114020 

USB JSB 20B,1 


0667 

01174 

177776 

MNUM OEC *2 


0666 

01175 

177776 

CNTON DEC -2 


0669 

01176 

OOOOOO 

NUMflR OEC 0 


0670 

01177 

037777 

MAX02 OCT 37777 

32768/2-1 

0671 

01200 

000100 

ONHUN OCT 100 

0672** 

****** 

*******1 

********** **********4 


0673** 

0674* 

PBE01CT09 algorithm CONSTANTS 

TAU*T/2».002HILLISEC0ND8 

0675* 

CHANNEL A 



0676* 

0677 

01201 

074000 

KlPAl DEC 30720 

»1«875(2**14) 

0678 

01202 

130000 

K2PA1 DEC -20480 

••1.25(2**14) 

0679 

01203 

014000 

K3PA1 OEC 6144 

■.S75(2**14) 

0680** 


*******4 
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FLUTTER SUPPRESSION SYSTEM COMPUTER PROGRAM (CONTINUED) 


0661* 

CHANNEL A FILTER CONSTANTS 


0662* 

0663* 

CONSTANTS FOR .002 

CYCLE TIME 


0684 

01204 

076566 

KlAll 

DEC 32118 

B«9e01960196(2**15) 

0685 

0686* 

01205 

000505 

K2A11 

DEC 325 

■,0099009901 (2**15) 

0667 

01206 

076566 

K1A12 

DEC 32118 

■.9S01960198(2**15) 

0688 

0689* 

01207 

000505 

K2A12 

DEC 325 

■.0099009901 (2**15) 

0690 

01210 

057234 

K1 A21 

DEC 24220 

*.7391304348(2**15) 

0691 

0692* 

01211 

010262 

K2A21 

DEC 4274 

■,1304347826(2**15) 

0693 

01212 

057234 

K1 A22 

DEC 24220 

*.7391304348(2**15) 

0694 

0695* 

01213 

010262 

K2A22 

DEC 4274 

*.1304347826(2**15) 

0696 

01214 

057234 

K1 A23 

DEC 24220 

*.7391304348(2**15) 

0697 

0698* 

01215 

010262 

K2A23 

DEC 4274 

*.1304347826(2**15) 

0699 

01216 

073076 

K1A31 

DEC 30270 

*1.847507236(2**14) 

0700 

01217 

106003 

K2A31 

DEC -29693 

=-.9061583110(2**15) 

0701 

01220 

021352 

K3A31 

DEC 8938 

*.2727596821 (2**15) 

0702 

01231 

001701 

K4A31 

DEC 961 

*,0293225116(2**15) 

0703 

01222 

160327 

K5A31 

DEC -7977 

■-.2434343704(2**15) 

0704*( 

»****«^ 


k*** ** *1 

Kr ********** kk k 


0705* 

0706* 

CHANNEL A OUTPUT SUMMATION CONSTANTS 

0707 

01223 

115657 

K0A12 

DEC -25681 

=-56.176/1.024*70(2**15) 

0708 

01224 

073350 

KOAll 

DEC 30440 

*66.508/1.024*70(2**15) 

0709 

01225 

152025 

K0A23 

DEC -11243 

*-24.593/1.024*70(2**15) 

0710 

01226 

027101 

K0A22 

DEC 11841 

*26,003/1.024*70(2**15) 

0711 

01227 

155422 

K0A21 

DEC -9454 

*-20.681/1 .024*70(2**15) 

0712 

01230 

077213 

K0A31 

DEC 32395 

*70.863/1,024*70(2**15) 

0713*^ 

k****«i 

k****«*4 

r** * * * *^ 

k***********« 


0714* 

0715* 

FILTER Initial conditions channel h 

0716 

01231 

000000 

Ylflll 

OCT 0 


0717 

0718* 

01232 

oooooo 

KK91 1 

OCT 0 


0719 

01233 

oooooo 

Y1012 

OCT 0 


0720 

0721* 

01234 

oooooo 

KXB12 

UCT 0 


0722 

01235 

oooooo 

Y1B21 

ncT 0 


0723 

0724* 

01236 

oooooo 

KY821 

OCT 0 


0725 

01237 

oooooo 

Y1622 

OCT 0 


0726 

0727* 

01240 

oooooo 

KX622 

OCT 0 


0728 

01241 

oooooo 

Y1023 

OCT 0 


0729 

0730* 

01242 

oooooo 

KXB23 

OCT 0 


0731 

01243 

oooooo 

Y1831 

OCT 0 


0732 

01244 

oooooo 

Y2831 

OCT 0 


0733 

01245 

oooooo 

X2B31 

OCT 0 


0734 

0735* 

01246 

oooooo 

X3B31 

OCT 0 


0736 

01247 

oooooo 

INPTB 

OCT 0 



B-17 



FLUTTER SUPPRESSION SYSTEM COMPUTER PROGRAM (CONTINUED) 


0757 

01250 

000000 

INPBl 

OCT 

0 


0738 

01251 

000000 

INP62 

OCT 

0 


0739 

01252 

OOOUOO 

OUTS ( 

OCT 0 

1 


0740* 


k******i 

ft ******1 

k * * * * 



0741*^ 

0742* 

0743* 

0744* 

» PPEOICTOR 1 
CHANNEL B 

algorithm constants 

TAUsT/2«.002HILLISECOND8 

0745 

01253 

074000 

KlPBl 

DEC 

30720 

al, 875(2**14) 

0746 

01254 

130000 

K2PB1 

DEC 

-20480 

a-l. 25(2**14) 

0747 

01255 

014000 

K3PBI 

DEC 

6144 

a, 375(2**14) 

0746*1 




***** 



0744* 

CHANNEL 8 FILTER CONSTANTS 


0750* 

0751* 

constants for .002 

CYCLE TIME 


0752 

01256 

076566 

KlBll 

DEC 

32118 

3.9801980198(2**15) 

0753 

0754* 

01257 

000505 

K2S11 

DEC 

325 

a, 0099009901 (2**15) 

0755 

01260 

076566 

K1B12 

DEC 

32110 

3.4601460146(2**15) 

0756 

0757* 

01261 

000505 

K2B12 

DEC 

325 

3,0099009401 (2**15) 

0756 

01262 

057234 

K1821 

DEC 

24220 

8.7391304348(2**15) 

0754 

0760* 

01263 

010262 

K2B21 

DEC 

4274 

8.1304347826(2**15) 

0761 

01264 

057234 

K1B22 

DEC 

24220 

3.7391304348(2**15) 

0762 

0763* 

01265 

010262 

K2022 

DEC 

4274 

8,1304347826(2**15) 

0764 

01266 

057234 

K1823 

DEC 

24220 

8.7391304348(2**15) 

0765 

0766* 

01267 

010262 

K2B23 

DEC 

4274 

8.1304347826(2**15) 

0767 

01270 

073076 

K1031 

DEC 

30270 

81,847507236(2**14) 

0766 

01271 

106003 

K2831 

DEC 

-29693 

a-. 9061583110(2**15) 

0764 

01272 

021352 

K3831 

DEC 

8438 

a. 2727598821 (2**15) 

0770 

01273 

001701 

K4831 

DEC 

961 

8.0293255116(2**15) 

0771 

01274 

160327 

K5031 

DEC 

-7977 

3-. 2434343704(2**15) 

0772*1 

»*****« 

^******^ 

k**** **i 

k **** 

******** 


0773* 

0774* 

CHANNEL 8 OUTPUT SUMMATION CONSTANTS 

0775 

01275 

115657 

K0B12 

DEC 

-25661 

8-56,176/1.024*70(2**15) 

0776 

01276 

073350 

KOBll 

DEC 

30440 

=66.588/1,024*70(2**15) 

0777 

01277 

152025 

K0823 

DEC 

-11243 

8-24.593/1,024*70(2**15) 

0778 

01300 

027101 

K0B22 

DEC 

11841 

826.003/1.024*70(2**15) 

0774 

01301 

155422 

K0821 

DEC 

-9454 

8-20.681/1.024*70(2**15) 

0780 

01302 

077213 

K0B31 

DEC 

32395 

870.863/1.024*70(2**15) 

0761*1 



k* ***** 1 

k *** * 

******** 


0762* 

0763* 

0764* 

0785* 

0766 

0787 

0768 

0764 

0740 

S REGISTER FLASH ROUTINE 
INITIAL CONDITIONS 

01303 177777 AFF OCT -1 

01304 177777 BFF OCT -1 

01305 177777 CFF OCT -1 

01306 177777 ACF OCT -1 

01307 177777 BCF OCT -1 

constants 

0741 

01310 

000000 

ERCNT 

OCT 

0 


0742 

01311 

177777 

ERCLK 

OCT 

-1 
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FLUTTER SUPPRESSION SYSTEM COMPUTER PROGRAM (CONCLUDED) 


0795* 

079«*CONSTANT3 

0795* 

0796 013i2 000002 AFFER OCT 2 

0797 01313 000020 BFFER OCT 20 

0796 01314 000200 ACFER OCT 200 

0799 01315 002000 BCFEP OCT 2000 

0600 01316 020000 CFFER OCT 20000 

0601 01317 177603 FLCNT DEC -125 

0602 01320 177777 FLAG OCT -1 
0603*********************************** 

oaoa* CROSS-CHECK initial conditions 

0605* 

0606 01321 000000 CHAER OCT 0 

0607 01322 000000 CN0ER OCT 0 

0606*********************** THE *** END 

0609 ENO 
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HARDWARE TEST PROGRAM 


0001 ASMB,A«BrT»L 

0002* FILE NAME BJTEST CREATED B/31/7B BY JRMATTHEN 

OOOS*******************************************ft***************************« 

0004* FR06RAM TO TEST ADC AND OAC CHANNELS 0*7 

0005* SET BIT IS FOR CONTINUOUS INFUT AND OUTPUT OF DATA 

0000* OR CLEAR BIT 15 FOR DISCRETE INPUT AND OUTPUT OF DATA 

0007* IN S REGISTER RHEN HALT EOB OCCURRS AND PUSH RUN 

OOOB* 

0004 10000 ORG lOOOOB 

0010 10000 103100 CLF 0 TURN OFF INTERUPT SYSTEM 

0011 10001 06200A LOA INIT INITIALIZE OAC 

0012 10002 102612 OTA 12B 

0013 10003 103712 STC 12B,C 

0014 10004 170140 INIT OCT 170140 

0015 10005 102020 MAIN HLT BOB SET OR CLEAR BIT 15 AS ABOVE INSTUCTIONS 

0016 10006 102501 LIA IB 

0017 10007 002020 SSA 

0016 10010 026013 JMP COUT JUMP TO CONTINUOUS ROUTINE 

0019 10011 026030 JMP OOUT JUMP TO DISCRETE ROUTINE 

0020 10012 102077 HLT 776 

0021 ********************************************************************** 

0022* ROUTINE TO INPUT AND OUTPUT DATA CONTINUOUSLY FROM 

0023* SELECTED ADC AND DAC CHANNELS 

0024* 


0025* 

HALT , 

22B INDICATES 

TOP OF 

ROUTINE— ENTER ADC CHANNEL INTO 

0026* 

BITS 

0*2 AMD 

DAC CHANNEL INTO BITS 3-5 AND PUSH RUN 

0027* 






0028* 

SET BIT 0 TO 

STOP CLEAR TO 

CONTINUE 

0029* 






0030 

10013 

102022 

COUT 

HLT 226 

READY FOR CHANNEL INFO 

0031 

10014 

016071 

JSB 

GCHAN 

PROCESS CHANNEL INFO 

0032 

10015 

016056 

LOOP 

JSB ADC 

GET INPUT FROM ADC 

0033 

10016 

062027 

LOA 

INPUT 


0034 

10017 

102601 

OTA 

IB 

PUT INPUT INTO S REGISTER 

0035 

10020 

072026 

STA 

OUTPT 


0036 

10021 

016044 

JSB 

DAC 

OUTPUT INPUT DATA TO DAC 

0037 

10022 

102501 

LIA 

IB 


0038 

10023 

000010 

SLA 



0039 

10024 

026005 

JMP 

MAIN 


0040 

10025 

026015 

JMP 

LOOP 


0041* 






0042 

10026 

000000 

OUTPT 

OCT 0 


0043 

10027 

000000 

INPUT 

OCT 0 


0044*4 

»***«*< 

*******4 

^****** 




0045* ROUTINE TO OUTPUT DATA IN S REGISTER TO DAC AND INPUT DATA 
0046* FROM ADC USING CHANNELS IN BITS 0-2 FOR ADC AND 3-5 FOR DAC 
0047* INPUT IS IN A REGISTER 

0048* 

0049* halt 21B MEANS TOP OF ROUTINE ENTER CHANNEL INFO 
0050* 

0051* halt 40B ENTER OUTPUT DATA 
0052* 

0053* HALT 31B SET BIT 15 TO STOP CLEAR TO CONTINUE 
0054* 

0055 10030 102021 DOUT HLT 21B ENTER CHANNEL INFO 

0056 10031 016071 JSB GCHAN PROCESS CHANNEL INFO 
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HARDWARE TEST PROGRAM (CONTINUED) 


0057 

10032 

102040 

HLT 

40B 



ENTER DATA 

0058 

10033 

102501 

LIA 

IB 



GET DATA 

0059 

10034 

072026 

STA 

OUTPT 




0060 

10035 

016044 

JSB 

DAC 



OUTPUT DATA TO OAC 

0061 

10036 

016056 

JSB 

ADC 



INPUT DATA FROM ADC 

0062 

10037 

062027 

LDA 

INPUT 



PUT INPUT INTO A REGISTER 

0063 

10040 

102031 

HLT 

31B 



SET BIT 15 TO STOP CLEAR TO CONTINUE 

0064 

10041 

002020 

SSA 





0065 

10042 

026005 

JMP 

MAIN 




0066 

10043 

026030 

JMP 

DOUT 




0067*« 


»******i 


»*****« 

r**** 

**** 


0066* 

OAC OUTPUT ROUTINE 

USING 

SLOT 

12B 


0069* 








0070 

10044 

000000 

OAC NOP 




0071 

10045 

066026 

LO0 

OUTPT 



»B*«OUTPUT 

0072 

10046 

101044 

LSR 

4 



SHIFT RIGHT 4 AND CLEAR UPPER BITS 

0073 

10047 

062111 

LOA 

DACCH 



'A»«OAC CHANNEL IN BITS 12-15 

0074 

10050 

030001 

lOR 

IB 



'A»*OV€RLAT OF OUTPUT AND CHANNEL 

0075 

10051 

102312 

SF3 

12B 



LAST OUTPUT complete? 

0076 

10052 

026051 

JMP 

*-l 



NO-CHECK AGAIN 

0077 

10053 

102612 

OTA 

12B 



YES-OUTPUT DATA 

0076 

10054 

103712 

STC 

12B,C 




0079 

10055 

126044 

JMP 

DAC, I 



RETURN 

0080*< 




^ ****« i 

Ift*** 

** ** 


0061* 

40C INPUT ROUTINE USING SLOT 

IIB 


0062* 








0083 

10056 

000000 

AOC NOP 




0064 

10057 

062110 

LOA 

ADCCH 



'A'anUTPljT CHANNEL CONTROL WORD 

0085 

10060 

102611 

OTA 

1 IB 



SELECT CHANNEL 

0086 

10061 

103711 

STC 

11B,C 




0087 

10062 

102311 

3FS 

1 IB 



INPUT complete ? 

0088 

10063 

026062 

JMP 

*•1 



NO-CHECK AGAIN 

0089 

10064 

10251 1 

LIA 

118 



YES-'A**INPUT 

0090 

10065 

012070 

AND 

mask 



mask out channel id 

0091 

10066 

072027 

STA 

INPUT 




0092 

10067 

126056 

JMP 

ADC, I 



RETURN 

0093* 








0094 

10070 

177700 

MASK 

OCT 177700 



0095*< 

k****«i 

»***«« *1 



k***« 

**** 


0096* 

GETS i 

IDC AND 

OAC CHANNEL 

NUMBERS 

FROM S REGISTER 

0097* 








0098* 

ADC CHANNEL 

IN HITS 0»2 




0099* 








0100* 

DAC CHANNEL 

IN BITS 3-5 




0101* 








0102 

10071 

000000 

6CHAN NOP 




0103 

10072 

102501 

LIA 

IB 



GET channel INFO 

0104 

10073 

070001 

STA 

16 



PUT INTO *8* 

0105 

10074 

012105 

AND 

MASKl 



*A'*eiTS 0-2 

0106 

10075 

032106 

lOR 

MASK2 



'A'slOOOON NbADC channel NIIMRER 

0107 

10076 

072110 

STA 

ADCCH 



STORE AOC CONTROL WORD 

0108 

10077 

060001 

LDA 

IB 




0109 

10100 

012107 

AND 

MASK3 



*A*«BITS 3-5 

0110 

10101 

001727 

ALF 

,4LF 




0111 

10102 

001200 

RAL 




*A* HAS OAC CHANNEL IN BITS 12-15 

0112 

10103 

072111 

STA 

DACCH 



STORE OAC CRONTROL MORD 
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HARDWARE TEST PROGRAM (CONCLUDED) 


0113 10104 126071 

0114* 

0115 10105 000007 

0116 10106 100000 

0117 10107 000070 

0116 10110 000000 
0116 10111 000000 
0120************** 
0121 


JMP eCHANfl 

MASKl OCT 7 
MA8K2 OCT 100000 
MA3K3 OCT 70 
AOCCH OCT 0 
OACCH OCT 0 

END 


RETURN 
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ALPHABETIC LIST OF INSTRUCTIONS FOR HP2100 COMPUTER 


ABS 

ADA 

ADB 

ALF 

ALR 

ALS 

AND 

ARS 

ASC 

ASL 

ASR 

BLF 

BLR 

BLS 

BRS 

BSS 

CCA 

CCB 

CCE 

CLA 

CLB 

CLC 

CLE 

fLF 

CLO 

CMA 

CMB 

CME 

COM 

CPA 

CPB 

DEC 

DBF 

DEX 

DIV 

DLD 

DST 

ELA 

ELB 

END 

ENT 

ERA 

ERB 


Define absolute value 

Add to A 

Add to B 

Rotate A left 4 

Shift A left 1, clear sign 

Shift A left 1 

'^And'^ to A 

Shift A right 1, sign carry 
Generate ASCII characters 
Arithmetic long shift left 
Arithmetic long shift right 
Rotate B left 4 
Shift B left 1, clear sign 
Shift B left 1 

Shift B right 1, carry sign 

Reserve block of storage starting at symbol 

Clear and complement A (Ts) 

Clear and complement B (I's) 

Clear and complement E (set E = 1) 

Clear A 
Clear B 

Clear I/O control bit 
Clear E 
Clear I/O flag 
Clear overflow bit 
Complement A 
Complement B 
Complement E 

Reserve block of common storage 
Compare to A, skip if unequal 
Compare to B, skip if unequal 

Defines decimal constants 
Defines address 

Defines extended precision constants 
Divide 

Double load 

Double store 

Rotate E and A left 1 

Rotate E and B left 1 

Terminate program 

Entry point 

Rotate E and A right 1 

Rotate E and B right 1 
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ALPHABETIC LIST OF INSTRUCTIONS FOR HP2100 COMPUTER (CONTINUED) 


EQU 

Equate symbol 

EXT 

External reference 

FAD 

Floating add 

FDV 

Floating divide 

FMP 

Floating multiply 

FSB 

Floating subtract 

HED 

P rint heading at top of each page 

HLT 

Halt 

IFN 

When N appears in Control Statement, assemble 
ensuing instructions 

IFZ 

When Z appears in Control Statement, assemble 
ensuing instructions 

INA 

Increment A by 1 

INB 

Increment B by 1 

lOR 

Inclusive "or" to A 

ISZ 

Increment, then skip If zero 

JMP 

Jump 

JSB 

Jump to subroutine 

LDA 

Load into A 

LDB 

Load into B 

LIA 

Load into A from I/O channel 

LIB 

Load into B from I/O channel 

LSL 

Logical long shift left 

LSR 

Logical long shift right 

LST 

Resume list output (follows a UNL) 

MIA 

Merge ”or" into A from I/O channel 

MIB 

Merge *'or” mlo B from 1/0 channel 

MPY 

Multiply 

NAM 

Names relocatable program 

NOP 

No operation 

OCT 

Defines octal constant 

ORB 

Establish origin in base page 

,ORG 

Establish program origin 

ORR 

Reset program location counter 

OTA 

Output from A to I/O channel 

0TB 

Output from B to I/O channel 

RAL 

Rotate A left 1 

RAR 

Rotate A right 1 

RBL 

Rotate B left 1 

RBR 

Rotate B right 1 

REP 

Repeat next statement 

RRL 

Rotate A and B left 

RRR 

Rotate A and B right 
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ALPHABETIC LIST OF INSTRUCTIONS FOR HP2100 COMPUTER (CONCLUDED) 


RSS 

SEZ 

SFC 

SFS 

SKP 

SLA 

SLB 

SOC 

SOS 

SPG 

SSA 

SSB 

STA 

STB 

STC 

STF 

STO 

SUP 

SWP 

SZA 

SZB 

UNL 

UNS 

XI F 

XOR 


Reverse skip sense 
Skip if E - 0 

Skip 11 I/O flag ■ 0 (clear) 

Skip If I/O nag » 1 (set) 

Skip to top of next page 

Skip if LSB of A « 0 

Skip If LSB of B » 0 

Skip if overHow bit « 0 (clear) 

Skip if overHow bit - 1 (set) 

Space n lines 
Skip if sign A - 0 
Skip if sign B » 0 
Store A 
Store B 

Set I/O control bit 
Set I/O nag 
Set overnow bit 

Suppress list output of additional code lines 
Switch the (A) and (B) 

Skip if A = 0 
Skip 

Suppress list output 

Resume list output of additional code lines 
Terminate an IFN or IFZ group of instructions 
Exclusive '*or'^ to A 
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